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PREFACE 


This study was performed as part of Contract Number NAS 5-11602 
entitled "Study of a Low-Altitude Satellite Utilizing a Data Relay Satellite 
System" for-the National Aeronautics and Space-Adminis-tration, Goddard 

Space Flight Center. The purpose of this study is to investigate the techni- 

cal considerations associated with a low-altitude satellite (LAS) operating 
in conjunction with a data relay satellite system (DRSS). One of the major 
problems is interference in the communication link between the two satellites - 
by transmitted radiation reflecte-d. and distorted by the earth' s-s.urface. This 
problem, termed "multipath," is the subject of th is rep ort. 
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INTRODUCTION 


1 . 


The descriptive term "multipath" refers to a communication link in 
which there exists more than one signal propagation path between the signal 
source and the receiver. The phenomenon of interest in this study occurs 
in the transmission of electromagnetic energy between-a low -altitude satellite— 
(LAS) and a data relay satellite (DRS) at geostationary altitude. For this 
study, low altitudes are in the 100 to 1000 mile., range. 

There are basically two paths between the satellites (Figure 1-1): a 
direct path, and a path by reflection from the earth's surface. The total 
reflecting surface may be thought of as being composed of many smaller re- 
flecting surfaces, and hence, the total reflection path is a collection of many 
reflection paths. The signal received via the reflection path is delayed and 
distorted with respect to the received direct signal, resulting in an undesir- 
able interference. Further, under frequently occurring conditions, the mag- 
nitude of the received reflected signal may be larger with respect to the direct 
signal power. Later in this report it is shown that, for a range of geometrical 
conditions and-a transmission frequency of 140 MHz over the sea, the average 
reflected power incident on the receiving antenna will be only 2 dB less than 
the direct signal-power. Interference of this magnitude has rather profound 
consequences on a comm unication system operating in this environment. 

The principal objective of this study is to quantitatively characterize 
the multipath phenomenon. Of major-interest is tile power received via both 
the direct and earth reflection paths. And, since free space propagation can 
easily be computed, the reflection process becomes the main subject of study. 

The goal, therefore, is twofold: . 

1) To develop a means for estimating the magnitude and nature of 
the reflected power 

2) To present numerical estimates of the reflected signal charac- 
teristics in graphical and tabular form for reasonable ranges of 
the important parameters 
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In addition, the fading characteristics of the combined signal are discussed, 
and a section is devoted to methods of communication improvement in the 
multipath environment. The major parameters includes 

1) Geocentric LAS-DllS separation angle, <p 

2) LAS altitude, h 

3) Transmission frequency, f 

4) Surface electrical properties (sea or land) 

5) RMS surface slope (roughness factor), i) 

6) RMS surface variation, o 

74 Antenna polarizations 

The reflection process is time varying and statistically nonstationary 
because of the changing geometry due to orbital motion*. However, the-short- 
term statistics may.be considered, stationary, thus allowing tractable 
statistical analysis of the reflection process. For a rough earth, the surface 
may be considered composed of many small patches, each with its own 
geometry, resulting in a reflected- power which is time delayed and frequency 
shifted.. Thus, the multipath phenomenon may be characterized by the 
following aspects of the total reflected power: 

1 ) Magnitude 

2) Time response of a transmitted impulse 

3) Frequency spectrum of a transmitted CW signal 
Quantitative estimates of these effects are presented-in Section 3. 

1.1 BACKGROUND 

At the heart of earth reflection analysis is the analysis of reflection 
from a flat rough surface, which did-not begin until about 1950. The analysis 
can be divided according to the type of surface model used. One type consi- 
ders surfaces composed of randomly distributed deterministic shapes. An 
example of such a model may be found in Reference 1 where Twersky analyzed 
a surface of hemispherical bosses distributed on a plane. A second model 
treats the surface height as a random variable. This approach seems more 
reasonable and has dominated rougu surface analysis in recent years. A model 
based on a statistical description of the surface roughness is used in this study 
and is discussed in Subsection 2. 4. 

Statistic a 11/, rough surfaces can be divided into three classes: 
slightly rough, intermediate, and very rough. ‘These classes can 
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bo quantitatively distinguished by means of a smoothness factor which 
includes the ratio of the rms surface variation and electromagnetic wave- 
length (see Subsection 2.2). But only the slightly rough and very rough 
surfaces can be treated mathematically with useful results; however, 
extrapolation between these two classes can yield useful estimates of sur- 
faces with intermediate roughness. Suppose that a perfectly smooth surface 
becomes slightly rough. The specularly reflected energy decreases as the 
roughness grows, but still contains the major portion of the energy. But, 
in addition, thcre-is a smaller amount of energy scattered away from the 
specular direction due to the roughness. The signal reflected specularly is 
phase -coherent, while the scattered signal is incoherent (random phase) at 
any distant reception point. As roughness increases, the coherent power 
decreases and incoherent power increases until, even in the specular reflec- 
tion direction, all received power is incoherent. The coherent component 
is computed using the expression for reflection from a smooth plane, with 
additional-factors included for the effect of roughness and the sphericity of 
the earth. The derivation of the expression is given in Reference 2, and the 
res ult-s -a . r - o .- d i s e u s s e d in Subsection 2. 3. 

Two approaches to the slightly rough earth are well documented. 

A perturbation technique— is used by Peake (Reference- 3) and Barrick 
and Peake (Reference 4) and later extended by Valenzuela (Reference 5). The 
results of these authors' approaches require the assumption that the rms 
surface variation be much smaller than a wavelength and that the Surface 
slopes be small. Another derivation technique employing the tangent plane 
approximation, or Kirchoff approximation, as it is sometimes called, is given 
by Davies (Reference 6) and Beckmann (Inference 2). These derivations 
require, in addition to the above a. sumptions, that the surface radii of 
curvature be everywhere much greater than a wavelength. The results of 
these analyses are useful for estimating the power reflected in directions 
other than the specular direction, but it can be shown that this incoherent 
power is significantly less than the coherent power reflected in the specular 
direction. 

Reflection from very rough surfaces has been analyzed using several 
techniques, all of which lead to equivalent results. Basic assumptions com- 
mon to these analyses are: 

1) The rms surface variation is larger than a wavelength 

2) The surface radii of curvature are much larger than a wavelength 

Assumption 2 is not always satisfied, particularly for land, but for the sea, 
which is a better reflector, it is satisfied quite frequently. Isakovich 
(Reference 7) provides the first comprehensive vector treatment of rough 
surface reflection from a physical optics point of view, and Davies essentially 
duplicated this work in a scalar formulation. This work was extended by 
Beckmann (Reference 2), Semenov (Reference 8), Stogryn (Reference 9), and 
others. Muhleman (Reference 10) considered the surface to be composed of 
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small, connected, perfectly conducting facets, and used the probability 
distribution of the facet normals in his analysis. Hagfors (Reference ) 1) 
extended these results, showing that they are equivalent to those of the other 
authors mentioned above. A third, but related, approach by Kodis 
(Reference 12) related the reflected power in a given c 14 icction to the number 
of specular points, i.e., the number of points satisfying specular geometry 
with respect to transmitter and receiver, Barrlck (Reference 13) 
extended this work by developing expressions for the number of specular 
points and average radii of curvature, obtaining results identical to 
the other analysis approaches. The form of these results Is discussed 
in Subsection 2, 4. 

The LAS/DLIS multipath problem has been considered by several 
investigators.. Durrani and Staras (Reference 14) began their analysis where 

the rough surface analysis discussed above ended. They took the general 

result for rough earth reflection,- which pertains to a surface patch, and 
formulated-the integral for total average received -power, The method of 
steepest descent was- used, along with a number of approximations, to arrive 
at a formula for the ratio of average received reflected power and received 
direct power which is defined as the "relative power" in Subsection 2,4. Some 
of the approximations made by Durrani and Staras are not accurate for 
all values of the geometric parameters. Further, because of the integral 
approximation approach,- an accurate assessment of the effects of antenna 
gain and polarization cannot be made... Massey (Reference 15) and Birch 
(Reference 16) are concerned principally with modulation, coding, and signal 
design to reduce multipath effects, and they add no additional information 
concerning the magnitude and nature of the reflected power beyond that given 
in Reference 14. 

In this study, the general rough surface solution of References 2, 8, 9, 
11, and 13 is used to compute the received incoherent power, and the 
slightly rough, coherent power is determined from the results of Beckmann 
a nd Spizzichino (Reference 2). The rough earth results arc achieved by 
numerical, integration of the surface integral via machine computation. This 
method numerically sums the power contributions due to many surface 
patches and, in the process, allows determination of the impulse response 
and frequency response of the reflection process. Thus, the multipath 
phenomenon may be characterized as mentioned previously. Further, the 
significant parameters may be readily varied and, in particular, the effect 
of antenna gain and polarization may be evaluated. 


1 . 2 PRINCIPAL RESULTS 

Estimates of the reflected power and the effects of the parameters 
listed previously are illustrated graphically and discussed in Section 3. 
Section 4 presents techniques for improving communication in the multipath 
environment. A few of the more significant results and conclusions arc 
presented here. 
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1,2.1 Interference* Magnitude 


In the absence of techniques for discriminating against the reflected 
signal, the interference by this signal can seriously degrade the performance 
of an LAB/DRSB communication link. If the electromagnetic radiation is 
polarized horizontally with respect to the reflecting area, i.e,, normal to 
the plane containing the LAS, DUB, and earth's center, then at 140 MHz the 
average magnitude of the reflected signal will be only 2 dll less than that of 
the directly received signal (Up -2 dD). This can be seen from Fig - 
ure 3-1. Further, from Figure 3-6 for this ease, the reflected signal power 
will exceed that of the direct signal approximately 20 percent of the time, 
and from Figure 3-7, the total received power will fade below the direct 
signal by 3 dD 10 percent of the time. 


1.2.2 Improvement Techniques 

A directional antenna on the LAS will provide a significant improve- 
ment by discriminating geometrically against the reflected power. However, 
many missions require an omnidirectional radiation pattern. For this case, 
circularly polarized antennas on both the LAS-suid DR8 will provide discrim- 
ination against the reflected signal. If an LAS antenna system can be designed 
to transmit/ receive the same sense of circular polarization in all directions, 
then the average reflected power will be less than 9 dB below the directly 
received power. This conclusion derives from Figures 4-3 through 4-6. 

Then, from Figure 3-7, the reflected power will be less than the direct 
power 99. 96 percent of the time and more than 3 dB below the direct signal 

98. 2 percent of the time. 

Signal processing techniques with improvement potential are dis- 
cussed briefly in Section 4, and include pseudo-noise (PN) coding, diversity 
techniques, frequency hopping, burst transmission, and data rate limiting. 

The most attractive and widely applicable of these techniques is PN coding. 
With this technique, the data signal is modulated by a signal corresponding 
to a PN sequence prior to transmission. This combined signal is demodu- 
lated at the receiver in a correlation process which discriminates against 
the reflected signal. It appears that effective use of PN coding requires a 
minimum RF bandwidth of approximately 20 kHz. Larger bandwidths will 
allow longer total communication time, better multipath rejection, and 
shorter signal acquisition times. 
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ANALYSIS AND MATHEMATICAL MODELING 


The multipath phenomenon requires mainly (geometric analysis in 
arriving at-a lathcmatical model useful for quantitative estimates of the 
reflected power. The analytical -treatment may be broken into throe funda- 
mental areas: 1) large-scale geometry involving the orbital positions of the 

two vehicles and the earth's mean surface, 2) smalle-r scale, surface rough- 
ness considerations, and 3) antenna radiation characteristics. These topics 
are treated in the above order in the following six subsections. 


2.1 LARGE-SCALE GEOMETRV 

In the following analysis, the LAS is considered as too uansmitter 
and the DRS as the receiver. This allows consistency and provides conven- 
ience in the ensuing discussions since the communication -link, including the 
multipath reflection effects, is bilateral. 

2. 1. 1 Separation An gle and Transmission Plane 

The two line segments from the center of the earth to the two- satel- 
lites may be used as a geometrical basis for analysis of the multipath 
phenomenon. The primary geometrical parameters are the length of these 
two lines and the angle between them, defined here as the separation angle, 

<p. In general, these three quantities vary with time unless the orbits are 
circular, in which case only the separation angle is varying. At any instant 
in time, the two lines define a plane which intersects the earth's surface. 
This basic geometry is illustrated in Figure 2-1. 

The direct transmission path lies in this plane, and the reflection 
from the earth's surface is nominally symmetric about it. This plane is 
subsequently referred to as the transmission plane, and it is a useful basis 
for the geometric analysis required in studying the reflection process. The 
separation angle may be related to the orbital parameters of the two satel- 
lites by considering the lines from the earth's center to the satellites as 
vectors in an inertial coordinate system, as shown in Figure 2-2. Referring 
to this figure, 
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Figure 2-1. Basic Geometric Quantities in 
LAS/DRSS Multipath Analysis 
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Once the time variation of 0, ^ , and X. is known, the time variation of <p may 
be determined from Equation 3. 

2.1.2 Coordinate Systems 

Using the transmission plane as a geometric basis, several useful 
coordinate systems may be defined. Referring to Figure 2-3, the coordinate 
system El E2 E3 is defined as follows: E] coincides with Rp), the line seg- 
ment to the DRS; E2 is orthogonal to El and lies in the transmission plane; 
and E3 completes the system, lying perpendicular to the transmission plane. 
All- vectors in the analysis of the reflection and direct transmission paths 
will be expressed in terms of this coordinate system. 

Two other coordinate systems required for describing the radiation 
patterns of the antennas are shown in Figure 2-3. The D] 1)2 1)3 system is 
parallel to the Ej E2 E3 system, but is centered at the DRS antenna. The 
Li L2 L3 coordinate system is centered at the LAS antenna. The Iq coordi- 
nate is parallel to the Rl line and points away from the earth, L2 is ortho- 
gonal to Li and is in the Ej E2 plane, and L3 is parallel to E3 and D3, 

2.1.3 Specular Point 

The specular point is that point on the earth's surface, in the trans- 
mission plane, where the angle of incidence is equal to the angle of reflection 
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measured from the local vertical. The geometry associated with this point 
is shown in Figure 2-4. It can be shown that 
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The specular point is the point where Gp, = e D» and is specified by the angle 
6 , which can be found as a function o Up by equating the two expressions of 
Equations 4-and 5. Figure 2-5 shows -6 as a function of <p for three different 
values of the LAS altitude. Figure 2-6 illustrates the relationship between 
the incidence angle, 0 , and the separation angle, <p . 

The specular point is at or near the center of the region on the earth's 
surface from. which the reflection toward the DRS occurs. Consequently, in 
summing the effects of reflection from many small surface patches, this 
point-will provide an origin for variation of the surface patch location. 

2.1.4 Scattering Patch Geometry and Reflection Plane 

Consider, a small patch, dS, of the earth's surface near the specular 
point- The location of the center of such a point can be specified by the two 
angles oi and (3, as shown in Figure 2-7- The angle o; is measured in the 
transmission plane, and (3 is measured in the plane defined by R a and E 3 . 
Thus, the area, dS, of the small patch is given by 


dS * R, 


cos 3 do; dp 


( 6 ) 


The vector Vl is directed from the LAS to the patch, and Vg) is directed 
from the patch to the DRS. These two vectors are given by 
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Figure 2-5. Specular Angles Versus 
Separation Angle 
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Figure 2-6. Specular Point Incidence 
Angle Versus Separation Angle 
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Also of interest is the vector, V B , which bisects the anple between V L and 

Vl> 
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The anple, Y, between Vg and Rg is piven by 
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The reflection plane is the plane containinp both Vj, and V j> Lnergy reflected 
from dS to the DRS propapates in this plane. 

2.1.5 Distribution of th e Sep ar ation Anple, y> 

Since the separation anple, <p , is the primary geometric parameter 
and several of the multipath characteristics will be determined as a function 
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of this angle, it in of into rout o * interpreting the data which is a function 
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Equation 3 relates <p to the of^tintical distribution of <p, some as sump " 

muflt be made * 
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For the purpose of numerical evaluation 
are made: 


the following assumptions 


1 ) 


between 0 and 360 degrees. 


» »»«fc*jwsar“ 


3) 


The inclination, i, remains constant. 


4) J [1C — 

The-last assumption determines ‘ from Equation 3 

be shown to be 1 08. 7 degrees. Any the earth blocks 

which is greater than this value is rr ^mngl separation angles. And 

transmission between the wo signal is received; 

when transmission 18 ^^inates or is reestablished via another DRS for 

fbnfi communication ter nunaves ox 1B . . 

which the separation unpin allows transmission. 


The altitude of the LAS is 500 statute miles. 


lie: r- 

Based on these assumptions, the Reparation angle, <P , 

relative frequency and the cumulativ { <\ The results are shown in 

for LAS inclinations of 0. 3°, 60. u £ an bc seen t h a t the density func- 

Figures 2-8 and 2-9. i rom g r r i c ^ re 2-9 can be used to estimate the 

tions peak at the inclination angle. gu^^ ^ lie8 in a given range. For 
percentage of time for w ic degrees, the separation angle, (p , lies 

example, for an inclination of-60 deg - _ * percent of the time, h igure 
between 40 and 80 degrees 66 .4 - X V 

2-10 indicates the effect of LAS altitude. 


2. 2 EARTH REFLECTION AND SCATTERING 


Mathematical modeling ^"n " 
earth's surface requires asbump - obtaining numerical results, 

arrive at analytical expre s sions j t more roughness conditions, 

Based upon the results ^corrospondinB £ ™° r °* con ditions may be esti- 

reflection characteristics foi ™ and roup h models are employed to esti- 
mated. In this study, botb 8 .... nd intermediate conditions. The 

mate reflection for these co . models will be discussed in the next 

with the nnnly.it. n.sminled with 
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Figure 2-8. Relative; Frequency 
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Figure 2-9. Cumulative 
Frequency of Separation 
Angle 
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thopo two extreme surface conditions, a discussion of what constitutes a 
rough or pmooth earth ip warranted, 


A. purfaee ip defined as smooth when the fields reflected by the flur- 
face from an incident plane wave propagate only in a single direction. v s 
i* iirto tint the fields reflected by one portion of the surface are highly cc r- 
Sod m InSxo rolaSoSp to fluid. .uattured by other portion, of the 
auS «o that the rofloctod onorgy propagate a. a wav. In a .ingle dtroc- 
tlon Similarly _a surface is qualitatively defined as rough when the fiol Is 
scattered by the' surface from an incidence plane wave are diffuse or propa- 
gate in various directions. 


To nr o vide aaauantitative criterion for surface roughness, consider 
two ooiat. on a flat surface separated by a height difference, h as shown in 
Figure 2-11. The waves reflected from these two surface points will diffei 
in-phase by an amount given by 


Al(j r.: 


4 it h COS 0 

TZ 


(13) 


where 0 is the angle of incidence as shown in Figure 2-1 1 and X is the EM 
wavelength. For a surface with randomly distributed surface variations, 
surface smoothness factor, q, may be defined by 


4 tt cr cos 6 
\ 


(14) 


where o is the standard deviation of the surface variation. Although the te is 
no sharp change from a smooth to rough surface, a criterion often u ®f cl a ® 
dividing^ line is the Rayleigh criterion which states that a mooth earth sat s- 
fies the following inequality 


IT 

2 


(15) 


If the surface variations are normally distributed, then it can be 
shown (see Reference 2, pages 80-89) that for a flat, perfectly reflecting 
surface, the fraction of the incident power which is reflected as coherent 
radiation is given by a coherency factor, p , defined by 


-q 


( 16 ) 
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Thus, a surface satisfying the Rayleigh criterion for smoothness may reflect 
only 8„ 5 percent of the ineUlont radiation coherently# 


A more restrictive smoothness criterion might require that the 
portion of the incident power coherently reflected by a smooth surface he . 
least 90 percent. This leads to-A-critorion for smoothness of; 


q < 0. 124 


(17) 


Similarly, a restrictive criterion for a rough surface might require that the 
portion of the incident power coherently reflected by a rough surface be no 
more than 10 percent. This leads to a criterion for roughness of: 


1 . 52 


(18) 


The region 0. 324 < q < 1.52 can be considered a transition region, both 
specular and diffuse scattering being present. 

The Beaufort sea state scale is presented in Table 2-1 where wind 
speed and the rms surface variation, <r,_are related. Figure 2- 1 2 _ shows the 
surface smoothness factor, q, as a function of <r and incidence angle for the 
VHF frequency of 140 MHz. The values of rms surface variation used for 
this figure correspond to those of Table 2-1. It can be seen that, for 
sea states the surface is not smooth except at very large incidence angle , 
and for most of the higher sea states, the surface is rough for a wrde range 
of incidence angle. But, for values of 0 between 0.03 and 0. 5 mrAcTB lhe 
surface roughness lies in the transition region, causing the leflected radia- 
tion, to have both coherent and incoherent components. !• or higher Irequcn- 
cies, the earth appears rougher and the reflected power will be principally 
incoherent for all sea states and all angles except near grazing. 

Modeling of rough surfaces based upon assumed surface statistical 
properties has been successfully atte mpted for three roughness conditions; 


1) q 

2 ) 0 


< 1 


3) q > 2. 1 4 

Barrick and Peake (Reference 4) develop expressions for reflected power cor- 
responding to 0 < q < 1 without using the Kirchoff approximation required by 
BecCann" (Reference 2) and others. However, as is noted in Re^n-4 
these expressions are most applicable to backScattermg, whereas LAS /DR o 
geometry leads principally to forward scattering except at small separation 
ancles. And. for the frequencies of interest in tins study, i.e. , f - 1 3(> MHz, 
the earth's surface will appear to be very rough (q ’ 2) most ot the time for 
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1 



Figure 2-12. Smoothness Factor, qj. 
for f * 140 MHz 


2-14 




" „tfL "JUW.- ’ -- 

y/ m 








TABLE 2-1. RELATIONSHIP AMONG BEAUFORT SEA STATE SCALE, 
WIND SPEED, AND RMS SURFACE VARIATION 


Beaufort 
Sea State 
Seale • 

Wind Speed, 
knots 

RMS Surface 
Variation, 
tr, .motors 

1 

1-3 

0. 02 

2 

4-6 

0.09 

3 

7-10 

0. 30 

4 

11-16 

0. 88 

5 

17-21— 

2. 1 

6 

2-2-27 

4. 0 

7 

28-33 

7. 0 

8 

34-40 

11.3 

9 

41-47 

17.4 

10 

48-55 

25. 2 

1 1 

56-63 

35. 0 


small separation angles. Thus, only two models are used in this study: the 
slightly rough model of Reference 2, which includes perfectly smooth condi- 
tions, and a wry rough model, developed and discussed in References 2, 8, 9, 
and 1*3. It is felt that intermediate conditions can be estimated from the 
results of these two models. Each model consists of a single mathematical 
- expression for the ratio of received reflected power and received direct 
power. For convenient reference, the models will be called -s-t nuni h and 
rough earth models, respectively. 


2.3 SMOOTH EARTH 

First consider reflection from a smooth plane surface. There is no 
Bimple rule for computing the exact fields reflected by such a surface, but 
an approximate field may be obtained as follows. When the Bource of the 
waves is far from the reflecting plane, the incident waves are substantially 
plane waves over any limited area. Because the angle at which a plane wave 
is reflected by the surface is equal to the angle of incidence, the reflected 
wave appears to come from an image source (Reference 2). Referring to 
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Figure 2-13. Reflection From a Smooth - 
Flat Earth 


+ = VERTICAL (IN PUNE) 
- = HORIZONTAL 



Figure 2-14. Polarization Vector 
Components Defined With ReBpect 
to Transmission Plane 


Figure 2" 13, the source or transmitter is located at T, radiating a field of 
strength S, The image is located at T' and of strength R- c times S, where 
R c is the surface reflection coefficient which depends on the angle of inci- 
dence and polarization of the wave. This is the same result obtained by ray 
theory or geometrical optics which arc based on the assumption that the 
electromagnetic fields boha-ve- locally as plane-siia-ves. 


Assuming for the moment that- the electromagnetic field is a scalar 
quantity, the equations for the— received direct and reflected signals are 
given, respectively, by 


P 


d s 


•> 


c L c 2 \“ 

2 2 
( 4 tt ) V c * 


P 


(19) 


°3°4 X 

(4 ,) 2 (V L + V D ) 2 


R 


( 20 ) 


where Gi is the gain of the transmitting antenna in the direction of the direct 
path; G2 is the gain of the receiving antenna in the direction of the direct 
path; P-is the transmitted power; X. is the wavelength; G3 is the gain of the 
transmitting antenna in the direction of the reflection path;.G4 is the gain of 
the receiving antenna in_the direction of the reflection path; R c is the Fresnel 
reflection coefficient; and Vl, Vpj, and V G are defined in Figure 2-13. 

Electromagnetic fields can be represented as vectors, and the reflec- 
tion of electromagnetic plane waves is dependent on the polarization of the 
EM wave where polarization refers to the orientation of the vector repre- 
senting the electric field, E. When considering reflection from a smooth 
surface, this direction is most conveniently referenced to the transmission 
plane, i. e. , the plane containing the incoming propagation direction, the 
surface normal, and reflection propagation direction. 

Referring to Figure 2-14, the polarization of an EM wave with respect 
to the reflecting surface and transmission plane is defined. The horizontally 
polarized plane wa.c'fi electric field is perpendicular to the transmission 
plane, while the vertically polarized component lies in this plane. Note, 
Figure 2-13 shows that, by convention, the component of E v i (the incident 
vertical field), which is parallel to the reflecting surface, changes direction 
upon reflection. All other field components remain in their original direc- 
tion. Using these conventions with Equations 16 or 17, only the Fresnel 
reflection coefficients are needed to compute the received power. 
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Figure 2-15. Magnitude of Sea 
Reflection Coefficients 


Figure 2-16, Phase of Sea 
Reflection Coefficients 








2.3.1 Reflection Coefficients 


A- theoretical expression for the reflection coefficient of a plane earth 
may bo obtained directly from Maxwell's equations. The earth is character- 
ized by its relative complex dielectric constant, c\ 


e' r -L- - 60 j\c 
c o 


( 21 ) 


where € is the earth's dielectric constant, c is the earth' s-conductivity, c 0 is 
the free space dielectric constant, and X is the free space-wavelength. Ihe— 
reflection coefficients of a sniooth plane earth surface are then given by: 

Vertical polarization: 


R 


v 


Horizontal polarization: 


R h 


g' cos 0 - 1 

fg' - sin^-G 

e V 



£* COS 0 + 

J e ' + sin^G 

COS 0 - ^ 

- sin^ 0 

COS 0 + . 1 

| 2 
g' + sin G 


(22a) 


(22b) 


where 9 is the angle of incidence. 

Figures 2-15 through 2-18 show the magnitude and phase of Rh and R v 
as a function of the incidence angle. For vertical polarization, the incidence 
angle where the rejection coeff? c ic n t - has - mln i mum m agnitude is called the 
Brewster angle. 

From these curves, the following points are of interest: 

1) The reflection coefficients for land are relatively independent of 
' radio wave frequency (for the frequency range of interest in this 

study). 

2) The pseudo-Brewster angle is 72 degrees for earth and between 
82 to 85 degrees for sea in the frequency range of interest. 
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Figure 2-17. Magnitude of Land 
Reflection Coefficients 
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Figure 2-18. Phase of Land 
Reflection Coefficients 




3) For angles of incidence less than the pseudo-Brewster angle, the 
angle of the horizontal reflection coefficient is near 180 degreos 
and the angle of the vortical reflection coefficient ip near 0 degrees. 
7 hu8 # for nnplcB of incidence than the pseudo -Prows ter 
angle, the horizontal and vortical reflection coefficients differ by 
approximately 180 degrees, which means a predominantly right 
circularly polarized wave will be reflected as a predominantly left 
circularly polarized wave. 

4) For angles of incidence less than the pseudo-Brewster angle, tho 

reflection coefficients of the so* are greater than 4hc -reflection 

coefficients of the land. 

2. 3. 2 Divergence Factor 

To this point, the model used in analyzing the reflection of EM waves- 
by the smooth. earth assumed the earth was an infinite smooth plane. A 
divergence factor, D, may be defined to account for the earth* s curvature*. 

It is defined as the ratio of. the power reflected by the spherical earth-to the 
power reflected -by an infinite plane. 

Referring to Figure 2-4, simple geometric calculations lead to a 
divergence factor given by (see Reference 2): 


( 1 + £ )( 1 + £/ cos 0) 


(23) 


where 


£ 


2v ls v ds 

R E < V LS + V DS ) 


(24) 


&£ is the radius of the earth. 

A more precise analysis by Bremmer (Reference 17) gives the diver- 
gence factor D as 


D 


1 

1 + i (1 + cos 2 G )/cos 0 + £, 2 


( 25 ) 


In later numerical computation, the more accurate Equation 25 is used. 
I-igure 2-19 shows the divergence factor as a function of separation angle 
and LAS altitude. h 
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2. 3, 3. Slightly Rough Surface 


It is shown in Reference 2 that, for a slightly rough surface where 
q « 1, a major portion of the incident power is reflected coherently in the 
speeular direction. The major additional assumption is that the radii of 
curvature of the surface be much larger than a wavelength, i. e. , that the 
surface consists of large-scale undulations. This condition is common for 
the sea, which presents the most serious multipath problems due to its 
highly reflective properties. If, in addition, the surface variation has a 
gaussian probability distribution, the power reflected by a slightly rough 
surface may be estimated by multiplying the smooth surface reflected power 
by a roughness coefficient which is i d e tv t i c a 1 - ■ ' to— the- c o h 0 r e n c y factor, p , 
defined by Equation 1 6 . 

2. 3. 4 Relative Power 

For evaluating the effect of earth reflection on the LAS/DRS. communi- 
cation link, the relative magnitudes of the direct and reflected power are of 
principal interest. The reflected power may be normalized to the direct 
power by defining the relative power, Rp, as the ratio of the two. Including 
the divergence factor and roughness coefficient, the re l at iv e .. . p owe r is given 
by 


P r , 

/ G 3 G 4\ 

1 

fM 

O 

i> 

i 

p d|r c | 

P d " 

V G 1 G 2/ 

+ V D>! 


Numerical results using Equations 16, 21, 22, 25, and 26 are presented 

in Section 3. 

Equation 26 represents the smooth earth model where slight roughness 
is included, i. e. , smooth here means not rough (q < 1 ). 


2.4 ROUGH EARTH 

T scattering of electromagnetic waves by a rough surface (q >> 1) is 
treated extensively in References 2, 4, 8, 9, 12, and 13. The specific prob- 
lem of interest in this study was also the subject of Reference 14. But, in 
the process of arriving at a formula for the relative power for the rough earth 
case, the authors made many approximations, some of which are not accu- 
rate for-all values of the geometric parameters. The approach in this study 
is to use the reflected power expression derived in Reference 2, but, whereas 
the authors of Reference 14 began their approximation of the required inte- 
gration starting with thiB formula, the intent here is to perform this integra- 
tion numerically via machine computation without approximation. 
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In addition to assuming that the smoothness factor, q, is large, a 
second important assumption in the derivation of the reflected power formula 
is that the radii of curvature of the surface are largo with respect to a wave- 
length. Although this condition is not always satisfied, especially for land, 
it is a common condition for the sea and for much land area. And since, as 
was mentioned in the previous subsection, the sea is the best earth surface 
reflector, this condition is not severely restrictive. 

To develop an expression for the expected value of the total received 

reflected power, P r , consider the power density, $ , ineiden t-ofl-a— small 

scattering patch, dS. $ is given by 


4> ■ 


C1 3 P 


4tt V. 


(27) 


where G 3 is the transmitter antenna gain in the direction of the scattering 
patch, P is the radiated power, and Vl is the distance from the LAS to the 
scattering patch, as shown in Figure-2.- 7. 

The expected value of the power reflected toward the DR S is defined in 
terms of a scattering cross section by 


<dP > = 4 > dC 
s 


(28) 


where dC is the scattering cross section for the scattering patch dS and is 
defined below. Using Equations 27 and 28, the expected value of the received 
power, Pj,, reflected from dS is given by 


<dP 

r 


g 3 g 4 \ 2 P 
_ _ - 

(4tt) J V L V D * 


dC 


(29) 


where G 4 is the receiver antenna gain in the direction of dS, and Vj) is the 
distance-from this patch to the DRS. As in the previous subsection, the quan- 
tity of major interest is the relative power, defined for a rough surface as 
the expected value (average) of the ratio of received reflected power to- 
received direct power. Combining Equations 19 and 29 
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whore S represent ft tho total pur fa c e a r o a common to tho view of both 
vehicles. 

2. 4. 1 Scatteri ng C ro ss-Section 

If the surface variation is considered to have random characteristics, 
whi ch appears to be a good assumption for a large area, then determination 
of the scattering cross section requires statistical analysis. The details of 
the required analysis are presented in Reference 2 and will not be repeated 
here. It is shown in Reference 2 that a statistical distribution of the surface 
variation about a mean value is sufficient for determining radiation fields. 

But this distribution does not describe tbe surface completely because it 
does not account for the separation of the peaks and valleys. A second func- 
tion, the autocorrelation function, may be used to describe this aspect of the 
Surface and is necessary for determining the average reflected power. In 
this study, the following two assumptions are made about the nature of these 
two function bj 

1) The surface height variation about a- zero mean value is distributed 
according to tbe gauss ian density function 


p( z ) 


1 _ 

(V \'2 TT 



where or is the standard deviation or rms variation value defined 
in Subsection 2. 2. For this study, this distribution is assumed to. 
apply over tbe entire surface of each elemental scattering patch. 

2) The surface autocorrelation function is given by 


A(x) 



where x represents the distance between two surface points, and 
T is a constant called the correlation distance. 


An important point concerning the autocorrelation function is that it applies to 
any two points regardless of the angular orientation of the line segment, 
between them; thus, the roughness is considered isotropic. The gaussian 
autocorrelation function is employed by Beckmann (Reference 2), Stogryn 
(Reference 9), Hagfors (Reference 11), Barrick (Reference 13), and others; 
but Fung and Moore (Reference IB) suggest a more complex function which 
they claim is better behaved. Although this function may be made to account 
for smaller scale roughness in addition to the large undulating roughness by 
introducing an additional parameter, it still behaves like a gaussian function 
for very large roughness q 1 . No experimental evidence has yet been 
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pro Rented to verify tbo autocorrelation function of Fung and Mooro, and, in 
SiSmon no Hinalo mathematical expression for the scattering cross. 

It Uon ia po^Mo. Thu™ tbo „J tractable ootutlon rculUnp from tbo 
above two statistical assumptions is used bore. 


Based upon this statistical description of the surface a lengthy dern a- 
tion, tbo details of which are given in Reference 2,. yields tbo following 
expression for tbo scattering cross section, dC: 


dC 



dS 


(31) 


where R~ is the reflection coefficient corresponding to 0 B in Figures 2-7 and 

^ » r .1 i.1. ^ 1 .x .. a 1 • r 1 1 /=» si 1 R ^ f) 6 Clll ftU/B 1 U 


2 -ZoTyTs the angle between Vjj; and tbc local vertical, R B , as shown in 


t IS me anuie » is* * . a r . i 

Figure 2-20; dS is tbe area of the surface patch; and \\ is defined as the 
roughness factor given by 


T 


The roughness factor can be shown to be the rms value of the surface 
slope with rospect to local horizontal. 


The scattering cross section dC corresponding to the surface patch dS 
can be interpreted as the portion of the patch dS which is oriented in such a 
way as toTrovide specula? reflection from the LAS to the DRS.. From Fig- 
ured 2-7 and 2-20, it can be seen that as the distance on the earth between 
dS and the specular point increases, the angle, Y , also increases and, hence, 
the scattering cross section decreases. 


2. 4. 2 Shadowing 


The derivation of Equation 31 assumes that every region oi the surface 
natch dS contributes to the scattered fields. This assumption neglects the 
shadowing of a surface by itself, an effect that is very important at large 
angles of incidence which correspond to large satellite separation angles. 
Indeed, the shadowing effect reduces the scattered fields to zero at an angle 
of incidence of 90 degrees. 


Several authors have treated the subject of shadowing. Beckmann 
(Reference 19) derives a shadowing function which yields the portion oi a sur 
face area not shaded, i.e., that portion which is illuminated by the incident 
radiation. This shadowing function involves the incidence angle, 0 , and 




■i fVypHi. 


2-27 


roughness parameter, i| , and it may be interpreted an the probability of a 
given point on the surface being illuminated. Beckmann claimed that this 
shadowing function should be included as a multiplicative factor in the 
expression for average reflected power given by Equation 29, 

However, Broekelman and Hagfors (Reference 20) claim that the 
shadowing function should be the conditional probability that a point is 
illuminated given that the surface at that point is properly oriented for 
reflection. Probability fundamentals favor this latter approach. The 
probability that the surface at a point is properly oriented is given by 
dC/dS|R c |“ (see Equation i I )» Thus, if 

A occurrence of properly oriented surface 

R. ■ occurrence of illumination (no shadowing) 

then the probability that the surface at a point will be both properly oriented 
and unshaded is given by 


P(A, B) P ( A ) P(B/A) 


where 


P(A) dC/dS|U c | 1 
P(B/A) shadowing function, S 

If the two occurrences were independent, then P(B/A) P(B) and 
Beckmann’s approach would be valid. 

Broekelman and Hagfors were unable to develop a shadowing, function 
appropriate to their approach, but both Smith (Reference 21) and Wagner 
(Reference 22) were successful. However, Smith considered only illumineu- 
tion of the surface by the transmitter, whereas Wagner included visibility 
by the receiver. If, 

C : occurrence of visibility by receiver 

Wagner's results yield the conditional probability P(D, C/A). Thus, for the 
LAS/DRS reflection signal, Wagner's results are more applicable. Wagner's 
analysis treats reflection in two dimensions, but the results appear to be 
applicable to the three-dimensional reflection process of LAS/DRS 
communication. 
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Several new quantities must be defined. The angles G^ and G^ are 
shown in Figure 2 = 20 and may be calculated from 


cos 0^ 
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( 32 ) 


(33) 


Now define 
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Two cases must now be considered in calculating the probability of visibility 
by both the transmitter and receiver. This probability may be expressed as 


P(B> C/A) = P(c/B, A) P(B/A) 


Case I: Forw ard Scatter 

Referring to Figure 2-20, if the angle between n^ and n r is less than 
90 degrees, then the reflection is called "forward scatter." lor this case, 
Wagner assumes that visibility by the receiver is independent of visibility by 
the transmitter. This assumption is consistent with that made by others in their 
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analyses. The expression for the shadowing function* S, is then calculated 
by Wagner to be given by 


u(cot 0.' - tan v' ) u(cot 0 r ' - tan y* ) 


t - 0 -‘ w ) 


(B. + D r ) 


where u( ) is the unit step function and 6 i' , 6 r ' , and V are the angles between 
V L , vd> Vg and the projection of R s onto the reflection plane. However* 
because of the geometric relationship which exists between Vl* Vd# R-S an< * 
Vg, tan di' > tan y* and tan 0 r > tan y and, hence. 
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1 - e 
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Case II:- -Ba.ck.:^.atter 

Referring to Figure 2-20, if the angle between n^ and nj. is greater — 
than 90 degrees, then the reflection is termed "backscatter" and 

P(C/B, A) - 1 


That is, if the point is visible from the LAS, then it is visible from the DRS. 
This is not strictly true for three-dimensional, reflection, but, for LAS/DRS 
geometry, this case arises principally when the incidence angles are small- 
So both P(C/B,A) and P(B/A) will be nearly equal to unity. The shadowing 
function then becomes 


u(cot 0.' - tan y ' ) \1 - e 1 
S - P(B/A) -- 1 B 


I 

where, as before, u(cot 6^ 


- tan y ) - 1 . 


Thus, 
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The criterion for distinguishing between the above two cases-is 

simple: 


n. • n 10 
1 r 


Case i 


n, • < 0 


Case II 



2.4.3 Divergence Factor 

It should be noted here that this rough surface analysis does not 
include a divergence factor. and rightfully so. The divergence factor is 
nroner for the smooth earth case, but for the rough earth, the earth s 
curvature is automatically accounted for by the ^e°m.tr»cal computations 
implicit in the determination of Y. V^.and V D in Equatio 


2.4.4 Relative Power 

Combining Equations 3 0 and 31 and including the shadowing func 
tion, the relative power-for a rough surface is given by 
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of R and the antenna gains 
c 


2. 5 DEPOLARIZATION 

Eauation 37 may be used to estimate the relative power for radia- 
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Figure 2-21. Shadowing Function 
S(0) = E(B,C/A) 
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Figure 2-22. Probah<Htie s of No Shadowing 
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perpendicular to the reflection plane, For instance, for vertical 
polarization, the reflection coefficient, R v , must be substituted for R c , and 
the gains Gj, G^, G3, and G, must correspond to the respective antenna 
gains for this, polarization, But, when the electric field vector is oriented 
arbitrarily with respect to the reflection plane, both reflection coefficients 
must be used in calculating the reflected wave, and the polarization of the 
reflected wave may be radically altered. This phenomenon is treated briefly 
in Reference 2, but the development and results are difficult to apply to the 
scattering problem of interest in this study. 

To derive the depolarization equations, the three coordinate systems 
shown, in Figure 2-23 will be used. The unit vectors are defined as follows: 
ki lies along the incident propagation direction; e hi is perpendicular to the 
plane containing k^ and the radius, Rg, from the earth's center to the reflect- 
ing point S; e v ^ is Orthogonal to k^ and. ej^j as shown; S3 coincides with the 
bisector of the angle between Vj and V r ; sj is perpendicular to the reflecting 
plane defined by and V r ; S2 lies in this plane c ‘■hogohal to sj and S3; k 
lies along the reflection propagation direction to the receiver; e hr is per* 
pendicular to the plane containing k^, and Rg; and lies in this plane per- 
pendicular to k r . These definitions may be translated mathematically as 

follows: 
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( 40 ) 


The electric field vector may be expressed in terms* of ' 
transmission unit vectors 


E. E. . e. . + E . e . 

i hi hi vi vi 


The quantities E } . and J£ y . are complex numbers defined by 


E, . Q, . e 
lu lu 


j(wt +<P h ) 


E vi = Q vi 0 


j(wt + </> v ) 


For- this analysis,- the time varying base, e , can be omitted since only 
the relative phases between the two components are important. Thus, 


E hi ■= Q hi e 


E . - Q . e 

VI VI 


The first step in the development is to decompose into components 
parallel to and perpendicular to s?. It is useful also to maintain the identity 
of the original polarization directions in these components. Thus, the fol- 
lowing definitions are made for the components of E. 


E hhs : P er P endicular to s 3* originated from E hi 


E hvs : P arallcl to s 3’ originated from 


E perpendicular to » 3 , originated from E yi 


E : parallel to s.., originated from E . 

V VS j Vl 
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It can be seen frdm the vector definitions of Figure 2-23 that 


E hhs = E hi [ (c hi * s i > ®] + (e hi *' s 2 ) s 2 1 

E hvs ’ E hi (e hi’ S 3 )S 3 

E , = E .1(6 . • s, ) s, + (e . • s-j) s 7 | 

vhs vi l ' vi 1 1 vi 2 2j 

E = E . (e . • s,) s, 
vvs vi vi 3 3 


( 42 ) 


In the reflection process, Efci ls and E v h s require the reflection coefficient 
for horizontal polarization, R h , while Ehvs and ^vvs require the vertical 
reflection coefficient, R v . The next step is to apply the reflection coefficient 
to the vector-field -components of Equations 42‘and then to express the field 
with respect to e^ r and-e v r coordinates. The- following definitions will 
s im plif y the resulting e xp r e s s ions . 



(43) 


Then 


( 44 ) 


E hhr = E hi (R h A hh + R v B hh» = E hi “hh 


E, = E.. (R. A, + R B. ) = E. . K, 
hvr hi h hv v hv hi . hv 


E . = E . (R, A . + R B . ) - E . K . 

vhr vi h vh v vh vi vh 


E = E . (E, A + R B ) = E . IC 
wv vi ' h vv v vv vi vv 


where E^, E v p R^, and R v are complex numbers. - The received electric 
field vector is given by 


E r = < E hhr + E vhr> e hr + < E vvr + E hvr> e vr 


= E, e, + E e 
hr hr vr vr 


(45) 


Equations 44 preserve the depolarization process in thatE^ vr is the compon- 
ent of the reflected electric field in the e v r direction , which originated from 
the e^i component of the transmitted field. 

2. 5. 1 Polarization Factor 

A convenient method for describing the field is via a polarization 
factor, P, defined by 


P = 



(46) 


since E v and Ej are components in some coordinate system, P is also refer- 
enced to that coordinate system. For instance, from Equation 41 


P. 

l 




j< W 


(47) 
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represents the polarization of the transmitted field with respect to the k^, 

e. .» e . coordinates. And 
hi vi 


E i = E hi (P i °vi + ' e hi ) 


( 48 ) 


From Equation 45, it can be seen that 


P r = E 


E + E, 
vr _ vvr _ hvr - 

hr E hhr + E vhr 


(49) 


Using Equations 44 and 46, P can be expressed as a function of P^. 


„ P . < R h A VV + R V V + R h A h.V + R V B hV F i K VV + V ..... 

r = P i < R h A vh + *v B vh> + R h Ahh + R v B hh ' P i K vh + “hh 


Since most of the reflected power is reflected from the earth's sur- 
face near the transmission plane defined in Subsection 2. 1, a special case 
of interest occurs for smooth earth geometry when V^, V r , and Rg of Fig- 
ure 2-23 are all coplanar, i. e. , when the incident and. reflected power both 
lie in the transmission plane. For this case e^i* Sp and e^ r are parallel 
and e v j, e vr , S3, and lie in the transmission-plane. It can be shown that 


B 


hh 


= A, . = B, ... = A , = B 


hv 


hv 


vh 


vh 


0 


A = -cos^ 6 
vv 


B = sin^ © 
vv 
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i ^ r 

ITT r*** ™ * 0 




'* 6 , 



and so Equation 50 reduces to 


P 


r 




0 


2 J 

COS 0 


(51) 


When 0 => 90 degrees, it can be seen from Figures 2-15 and 2-16 that 
R v » R. ; thus, P r <s P i# and the -two coordinate systems coincide. Gonse* 
quently, the polarization is unchanged. When 0*0 degrees, P r - -P^, 
thus the reflected polarization is 180 degrees different in phase with respect 
to the receiver coordinates than that which was transmitted. Then, for 
instance, if the LAS is transmitting a right circularly polarized wave omni- 
directionally away from the spacecraft (a-hypothetical case for illustration, 
only), and the DRS is "overhead" (<p = 0, 0 « 0), then the wave reflected 
toward the DRS will be left circularly polarized. This phenomenon allows 
the possibility of discrimination between the dixjo.cL.and- reflected signals and 
is discussed further in Section 4. 


2. 5. 2 Power Reflection 


For a general polarization, a power reflection coefficient. F r , can 
be defined. This reflection coefficient replaces the quantity JH-cT * n Equa- 
tions 26 and 37 and is defined by 1 


r 


r 



specularl y reflected power 
incident power 


E 


hr 


E hi 


E 


vi 


E . 
vi 


From the definitions of the polarization factors given by Equations 47 
and 49 


E hr 

2 

(1 + 

P r 

2 ) 

E hi 

T 

(1 + 

P. 

i 

2 > 


( 52 ) 
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6 ANTENNA CHARACTERISTICS AND RECEIVED POWER 

The characteristics of interest here are the radiation pattern and 
polarisation which must he related to the scattering process. The polarisa- 
tion factor described in the previous subsection will be useful in represent- 
ing the gain and polarisation of the antennas. There are five polarisation 
factors of interest which are referenced to coordinates shown in figure <£-<-.-• . 

P. * incident factor referenced to c hl » P yl # hj 
P - reflected factor referenced to e hj ., e^, k y 
P s transmit direct factor referenced to e^, e^, 

p s receiving antenna direct factor reference to e^, e v( j* 
rd 

P - receiving antenna reflection factor referenced to i> hr , o yv » -k r 
t V 

From the orbital geometry, it can be seen that vectors and k r are 
separated by a small angle. It is shown in Section 4 that this angle' is 
always less than 3 degrees for LAS altitudes less than 800 n. mi., Thus, i 
the DRS is the receiving satellite, P f( jw I 3 rr > and 18 tllc transmitter, 

P d -P £ . 


2. 6. 1 Received Power 

Considering only specular reflection and depolarization, the ratio, 
R q , of received reflected power to direct power is given by 


E hr E hrr * E vr E vrr 
L ° = E hd E hrd + E vd E vrd 


E. 2 E 2 
hr hrr 


E hd E hrd 


IIP P 
r r rr 

2 1 + P d P rd 2 


Now, by definition of the antenna power gains Gj, G 2 , G 3 , and G 4 in 
Equations 26 and 37, 


c i = E h.J 2 + l E vdl 2 " l E i.dl 2 (' + KP 


^'2 ' E hrd + E vrd E hrdl (* + P rdl 
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E 


hll 


E. 


vi 


E hi | 2 (' + p t | 2 ) 


G, 


E 


hrr 


E_ 


vrr 


E 


hrrl 


(' * p ,, 2 ) 


Using Equations 39 and 46, Equation 53 becomes 

/ C 4 KrIM 1 IfrZr/l! (j * l P dH ( l + l P rdP) 


R 


o y g 1 g 2 


1 +P d P rd*l“( 1 + l P «| 2 ) 


From Equations 44 and 45, 

E. = E. . K. . + E . K . - E, . (K. , + P. K .) 
hr hi hh vi vh hi ' hh i vli 

Then from Equations 54 and 56, 


G~ G. 

R -r ~~tA IK, , + P. K 


o “ Gj G 2 r'hh T A i vh 


PC 


where PC is the polarization coefficient defined by 


PC =* 


l + p r p tr 

*(>♦ 

p d 

I 2 ) (> ♦ 

P rd 

2 ) 

1 + p d p ,d 

*(>♦ 

p. 

1 

2^ 

i 

H 

(*♦ 

P rr 

2 ) 


( 54 ) 


(55) 


(56) 


(57) 


(58) 


For linearly polarized antennas, P^, P^, P r d* anc * ^rr are positive real 
numbers. Horizontal polarization corresponds to P ~ 0, while vertical 
polarization corresponds to P - ®, which presents computational problems 
in evaluating R 0 as given by Equation 57. To overcome this difficulty, the 
inverse polarization factor, Q, may be defined by 


Q - 


P 



(59) 


It can then be shown by a derivation identical to the one above that 


It 


' G 3 G 4 ’ 


o l Q,0. 


K + Q. K. 

vv l hv 


QC 


where 


1 +t V^rr 

KnJ 

ik 

z ) 

(i + 

Q rd 

2) 

1 + Q d Q rd 

1*0 ♦ 

Q i 

z ) 

6 + 

Q rr 

F) 


It may be noted from Equations 57 and 60 that 


QC = PC 


K. . + P. K . 
hh i vi%- 

K + Q. K. 
vv i hv 


( 60 ) 


(61 ) 


(62) 


Then for horizontal polarization, P - 0, and for vertical polarization, 
Q = 0. 

2. 6. 2 Relative Power 


In Equations 57 and 60, the terms [k^ + Pi K v ^| 2 and 
K vv + Qi Kj lv p are equivalent to | R c | 2 of Equations 26 and 37. The 
polarization coefficient PC or QC is an additional term that must be inserted. 
Thus, expressions for the relative power for both the smooth and rough earth 
cases can be made to include the antenna characteristics and the effects of 
depolarization. 

2. 6. 3 Rough Earth Case 


R 


_ v i_ f 

^ Js 


C 3 C 4 


V, 


D 


K, . + P. K , 
hh t vh 


PC 


2 4 

H cos y 


8(G) exp 




( 63 ) 


where dS is a small scattering patch on the earth's surface and each of the 
terms on the right side of the integral sign correspond to each dS. 


1 


2, 6, 4 Smooth Earth Cane 

Home simplifications can be made for the smooth earth ease by noting, 
an wan done in Subsection 2, 5, that K^i Hi,, K = R v nin^h - Rp con^O, 
and K hv - K vh - 0. Then 



Equations 63 and 64 are the basis for some of the relative power calculations 
discussed in this report. 

2. 6. 5 Antenna System Co ordinates 

The above development with the resulting expressions requires that 
the LAS and DRS antenna characteristics be known for every significant scat- 
tering point as well as for the direct path direction. For actual calculation, 
the gains and polarization factors can be expressed as functions of transmis- 
sion direction with respect to a Coordinate system fixed on each satellite. 
Such coordinate systems were defined in Subsection 2. 1 and are shown in 
Figure 2-24. 

The two angles 4 ^ and <p relate the transmission direction Vj^ to 
the LAS coordinate system. Referring to Figure 2-24, it can be shown that 





whore sgn( ) means "sign of" and the numeric subscript refers to the 
component in the Lj b z L .3 coordinates of Figure 2 = 24. And, since reflection 
is symmetric about the transmission plane, 2 d, and so 

-90 degrees £ <pj & 90 degrees 


0 < . £ 1 80 degrees 

L 


Similarly, for the DRS antenna coordinates: 


V - D 
V D 1 


cos <p. 


r D3I 

V D | smC D 


Mr v -W 1 D3I 

= s en(V D2 )cos f |V D | sinCD 


The gains and polarization factors can be expressed as a function of £t,' <t»L* 

£ d , and <p D after the antennas are related to the L 1 L z L 3 and Dj D 2 D 3 
coordinate systems. 

2.6.6 Horizonta l and Vertical Polarization 

To provide numerical results useful for. under standing the effects of 
the several parameters, two linear polarization orientations will be used. 

The horizontal polarization vector for both satellites will be taken as the 
normal to the transmission plane (the plane containing the LAS, DRS, and 
earth's center.). The vertical polarization direction for each satellite will 
be parallel to the transmission plane, but an arbitrary choice of direction is 
required. 

For the LAS, an omtii capability for vertical polarization will be 
assumed. That is, the gains in the direction of the specular point and in the 
direct path direction will be taken as equal. This will allow the reflection 
process to be characterized nearly independent of the antenna characteristics. 
Thus, for any given value of the separation angle, there are two equal vertical 
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polarization vectors, Vjj; ancl Vj^d* corresponding to specular point incidence 
and the direct path. These two vectors are illustrated in Figure 2-25 from 
which it can be shown that 


V 


LI 



( 69 ) 


where 


a 90 - ^ - (P 


and the E^ E^ coordinates of Sub-slcct ion 2. 1 are used. 

For the DRS, the vertical polarization direction is chosen as the 
perpendicular to its earth radius vector, as shown in Figure 2-25. Thus, 


V 


PD 


/ 0 \ 

1 

\ 0 / 


(70) 


This appears to be the most reasonable since the DRS antenna is most likely 
to be an earth coverage antenna pointed, nominally, at the center of the 
earth. Using these coordinates, the horizontal polarization vector, H, 
defined above can be expressed by 


H 


/ o\ 

0 

M / 


(71 ) 


With these definitions, and considering separately, horizontal and vertical 
transmitters with 1 watt of power, the gains and polarization factors could 
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T 



bo defined according to subsections 2, 6, 1 through 2,6,4, However, this 
case is of sufficient interest that additional relationships will be developed, 
First, define the following quantities : 


G tn , II • 0 , . ' c, . , 

Hhi hi hi.i 


G . II • e . ” e ... 
Hvi vi vii 


(72) 


G, rl . - V,-, T • e, . ~ e, . . sin a + o. . cos a 
Vhi PI, hi hil hi 2 


G. r . - V,., T * 0 .- 0 ., sin a + e ... cos a 
Vvi PL vi vil vi2 


G = H • e - e 
nHr “ e hr hr3 


G TT = H • e „ = e 
vHr vr vr3 


G hVr = V PD ‘ e hr “ e hr2 


G , r = V,.,_ • e - e 
vVr PD vr vr2 


where the numerical subscripts corresponds to vector components in the 
Ej E^ E^ coordinate system. 

It can now be shown that, if horizontal polarization is transmitted, the 
reflected, horizontally polarized field is given by 

E Tf „ TT ■= G. „ (G tlU . K, . + G tI . K . ) + G „ (G„. . K, + G tI . K ) (73) 

IIllH hill* ' Hhi hh Hvi vh vHr ' Hhi hv Hvi vv ' 


where K vv , K v j a , Khv» and K vv are defined in Equation 44. The reflected 
vertical polarization is given by 


E tI1) ,, G, ,, (G it . . K. , -!• G . . . K . ) + G .. (G.- . K. 4 G tI . K ) (74) 

HRV hVr Ilhi hh Hvi vh vVr Hhi hv Hvi vv 
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Similarly, if vortical polarization in transmitted, the received 
reflected vertical and horizontal polarization is given respectively by 


E 

f VRV 

E V11H 


Ci. .. (G ... , 
hVr Vlu 

G hllr (C, Vhi 


K, . + C.„ . K . ) + G TT (G. n . JC + G . K ) 

hh Vvi vh vllr Vhi nv Vvi 


vv 


K. , 1“ G. r . K , ) 4* G TT (G am , K, + G. r , K ) 

hh Vvi vh vllr Vlu hv Vvi vv 


(75) 

(76) 


For direct transmission, the received horizontal polarization is given 
by Equation 19 with G{ ~ G? 1. For vertical polarization with the con- 
vention chosen above, Gj 1 and G 7 -- cos^v, where v is defined in 
Figure 2-25. 

Combining these results with the scattering cross section as was 
done above, the relative power for transmitting and receiving horizontal 
polarization (horizontal-horizontal) is given by 


V 


(R 


P 


HH 


4tt 



2 4 

\\ cos y 


exp 


'-tan 2 Y > 


IdS 


(77) 


The horizontal-to -vertical relative power is given by this equation with 
^HRV substituted for 

The vertieal-to- vertical relative power is given by 


<V 

vv 



4tt cos v 


E 


VliV 


s 2 (0) 


2 2 2 4 e X P 

'S V_l V D n cos Y 



(78) 


and vertical -to -horizontal is given by substituting f° r ^yRy * n 

equation. 
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2,6,7 Smooth Earth 


m 


m 

m 

sa? 

'w 

=§ 




Mv#v 

-Vr*| 


For the smooth earth model where reflection comes from the 
specular point. 


Is 

> 

H. 

1 1K IT 


h 

1 

2 


VRV 


R 

V 


> » 
cok“t 


where t is defined in Figure 2-25. 
Thus, 


<V 


IIH 


(R P > vv 


V 


V -I- V 
V LS DS> 


D 


R, 


V 


Vs + Vs 



COS T 
C OS V , 


D 


R 


( 79 ) 


(80) 


(81) 
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3. MULTIPATH CHARACTERIZATION 


The multipath phenomenon cannot easily be described by one quantity 
or function. However, it can be usefully characterized by its major effects 
upon a transmitted signal. The previous section developed expressions for 
the relative power (ratio of received reflected power to received direct- 
power) for both the smooth and rough earth cases invoiving Jke™ a J or -- 
parameters of satellite-to- satellite communication and the reflection 
process. These expressions may be used to evaluate the following 
three distinct but: interrelated aspects of the multipath phenomenon. 


1 ) 


Power transmission - relative power reflected by the 
earth to the receiver 


2) Time response - rec 
time 


eived reflected power as a function of 


3) Frequency effects - multipath effect on a fixed frequency 
signal 

Ir. providing quantitative estimates of the above multipath aspects, the 
effects of the following major parameters will be shown: 


n Geocentric LAS-DRS separation 




2) LAS altitude, h 

3) Transmission frequency, f 

4) Surface electrical properties (sea or land) 

5) RMS surface slope (roughness factor), i) 

6) RMS surface variation, a 

7) Antenna polarizations 



SEPARATION ANGLE, DEGREES 
a) Horizontal Polarization 



b) Vertical Polarization 

Figure 3-1* Relative Power Cor Nominal Multipath Parameter 



r 


(roez-eisoo 



3. 1 POWER TRANSMISSION 

Power transmission has three aspects of interest: 1) magnitude 

of the received average reflected power with respect to the received dnect 
power, 2) statistical variation about this average which results in a fading 
signal power, and 3) distribution of reflected power over- the earth's 
surface. These three aspects are discussed in this subsection. 

3. 1. 1 Reflected Power 

A useful quantity for evaluating the received reflected power with 
respect to the received direct power is the relative power defined, in Sub- 
section 2.3. The relative power, Rp, is defined as the ratio of the. average 
received reflected power to the received direct power.. For the smooth 
earth model of Subsection 2.3, no statistical averaging is necessary, but for 
the rough earth model, such averaging is required as discussed in 
Subsection 2. 4. 

The final mathematical expressions are given by Equations 63 and 64. 
These equations require that the gains and polarizations of both LAS and 
DRS antennas be known for the direct path and over the scattering region. 

For purposes of characterizing the multipath phenomenon, two orthogonal 
linear polarizations were chosen and are described in Subscction-2. 6. 6 where 
the relative power expressions for these two polarizations-aro derived. 

Using these equations, a digital computer program was written to 
compute the relative power for the smooth and rough earth models as a 
function of the above mentioned parameters. The program and details 
associated with the numerical computation arc discussed in~Appendix A. 

Since the separation angle, <p, is the major geometric parameter, the fig- 
ures of this subsection, will plot relative power as a function of cp with the 
effects of other parameters shown on these coordinates. 

Surface Type 

It can be seen from Figures 2-15 and 2-17 that reflection coefficients 
for the sea are greater than those for land. The difference between the 
reflection coefficients is due to the different electrical properties of the 
two surfaces, which are accounted for by the relative dielectric constant 
€ j. = e/e Q a.nd conductivity, c, in Equation 21. For this study, the following 
two sets of parameters were chosen based on Reference 23 (Table 5. 1, 
page 398) which arc representative of each surface type: 

Sea: c . - 80 c = 6 mho /meter 

Land: € - 10 c 10 mho/meter 

For sea reflection, Figures 3- la and b show the relative power as a func- 
tion of separation angle, <p, for horizontal and vertical polarization, 
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RELATIVE POWER, DECIBELS 
















respectively, and a transmission frequency of 140 MU'? , Both the perfectly 
smooth earth (a s Q) and rough earth (q s 0. 1) case a are shown along with 
the coherency factor, p, defined in Subsection 2. 2. The coherency factor 
corresponds to rr = 1,0 meter. As mentioned in that subsection, the transi- 
tion region can be considered to correspond to a range in the cohe-reney 
factor, p , between 0.1 (“10 dB) and 0.9 (-0,457 dB). This region is shown 
in Figures 1-la and b. Note that the results for the perfectly smooth and 
rough earth cases are nearly identical, even in the transition region; thus, 
it is not unreasonable to assume that the relative power will lie between 
these two values in the transition- region or, at most, be equal to their sum. 

Equally important is the nature of the reflected power. As discussed 
in Subsections 2.3 and 2.4, the power reflected from a smooth or slightly 
rough earth is phase coherent, while that reflected from a rough earth is 
incoherent. In the transition region, it is to be expected that the total 
reflected power will consist, of both typos. Using the coherency factor cor- 
responding to the given value of o-, the relative incoherent and coherent 
power, respectively,., are shown in Figures 3-2a and b,. This method of 
representing the relative power appears to bo the most appropriate and 
will be used subsequently in this- -suction,. 

Figures 3-2a and b also show the relative incoherent and coherent 
power for nominal land reflection. Since the reflection coefficients for land 
are relatively independent for frequencies in the range of 140 MI! a to 
15 GHz, the incoherent curve remains the same over this range, but the 
coherent power is reduced to a negligible amount for nearly all separation 
angles with frequencies above 2 GHz. It cart be seen that the relative power 
for-sea-reflection is greater than that for land reflection. 

Inasmuch as the sea covers approximately .’0 percent of the earth's 
surface and, at the same time, represents the worst-case reflection, this 
surface will be used for the remainder of this section to show the effects 
of other parameters. Typical roughness parameters for the sea correspond 
to 0. 05 < q < 0. 1 and it* 1,0 meter. When the effects of other parameters 
are being studied, q = 0. 1 and <r - 1.0 will be used. 

Transmission Frequency 


The effect of transmission frequency can be seen in Figures 3-3a 
and b. As the frequency increases, the earth appears more and more 
rough, with the result that the reflected power is incoherent for almo-st all 
separation angles at frequencies above 2 GHz, 
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LAS Altitude 


The effect of LAS altitude on the relative power is shown in 
Figures 3 -4a and b for a frequency of 140 MIL.. The relative power 
decreases an the altitude increases, 

Roughness Parameters 


The roughness factor, q, representing the rms surface slope and the 
surface rms height variation, a, are clearly not independent variables as 
can be seen from the definition 


n 


z o' 

T 


where T is the correlation length. However, as the surface becomes rougher, 
it. is apparent that o- increases and T decreases. Thus, q will experience a 
greater variation than a. Experimental studies discussed in Reference 2 
show that the rms slope varies between 0 and 16 degrees, corresponding to 
a maximum variation in q from 0 to about 0.3. More common values cor- 
respond to 0.05 s q £ 0. 1. The following data indicate the relationship 
between q and ir used to produce Figures 3-5a and b. 


n 

cr, 

meter 

0. 05 

0. 05 

0 . 1 

1.0 

0. 2 

1. 5 

0. 3 

2. 0 


3. 1 . 2 F ading 


It is of interest to consider the fading of a single transmitted CW 
frequency. Modulation on this signal will change the fading characteristics, 
but the analysis of this situation is beyond the scope of this study. A 
coherently reflected signal will cause fading with a beat frequency equal to 
the difference in dopplcr shift between the direct and specular reflection 
paths. This difference depends on both the separation angle and satellite 
velocities with respect to the transmission plane. Doppler shift is dis- 
cussed in Subsection 3.2. If the reflected and direct signal powers are 
nearly equal (Kp « 0), then the fade will be deep. The ratio of total power 
and direct power will have the maximum and minimum values given by 


fMP 2KP >* 


figure 3-6. 

o£ Total Received Power 



_ Ornbabi' ^.y That Receive 
.^fuGr'eater Than K Time* Average 


d Reflected 

Power, <1 


( 0 ) 815 -^ 0 ®° 




P, ! P 

(I l* 


n 


i> 


1'hls relationship in shown in Figure i-6, Thu coherent case is determinis- 
tif and easily calculated tor a particular geometry and a value for K . How- 
ever, l lu* incoherent cany requires considerable analysis. 


In Subsection 2.4, the relative power tor incohc r ciil reflevdion from 
a rough earth surface is defined as the ratio of the average reflected power 
and direct power. Variation from this average results in fading. Un de r the 
assumptions about the statistical nature of the surface given in Sub- 
section 2.4, it can be shown that the envelope of the reflected electric, field 
has what is known as the Rayleigh distribution (Reference 2). The voltage in 
the receiver is directly proportional to the- ole civic field, and so the prob- 
ability of the envelope of the voltage, v, lying between I lie- values a and b is 
given by 



dx 


where r is the rms value ol the voltage. Since the received power is 
proportional to the square of the voltage, a change of variable yields the 
distribution for the received reflected power P t . 


P(P r 2 A) 



:P > 
r 


dz 





where < P v ^ denotes the average reflected power (see Equation 29). Thus, 
the probability that the instantaneous rutlooted power is K times greater 
than the average is simply given by. 


P(P r ? k <P '•) 


and is plotted in Figure 3-7, Note that the probability of P r being greater 
than v.p |t > is only 0.37, not 0. 4, 
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Figure 3 


X, DECIBELS 

-8. Probability Distribution of Ratio, Pjj_, of Total 
Instantaneous Power and Direct Power 



PROBABILITY THAT ORDINATE VALUE WILL BE EXCEEDED 


Figure 3-9. Inverse Probability 
Distribution of P» 
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If the total — reflected plus direct signal '-is considered, the 
instantaneous power, Pj., can be shown to have the Rician-disiribution 
given by 


-P, \ 2\ r P 


p(P t ) = exp \-^p 


d \ t 

> I <p > 

v / r 


• P. 


•> C>X P l<p > 


Defining the power ratio, P^, by 


/^ p t p d 

\ <P r" 


P 


R 



where, again, P fc is the instantaneous total power. The probability distri- 
bution of Pr can c be determined as was done in Reference 24. The results 
are repeated here as Figure 3-8. 


This figure may be used to estimate the amount of time in which a 
given. amount of fade will occur. For instance, if Rp - -4 dB-, then from 
Figure 3-8, the amount of time that P^ < P^, i. e. , P(Pr < 0 dB-), is about 
40 percent. And the probability that P(; < P^/IO is about 2 percent. The 
inverse probability distribution-is shown in Figure 3-9. A special proba- 
bility scale is used to show the probability that Pj^ will be greater than a 
given value. From this figure for R„ - -4_dB, it can be seen that 
95 percent of the time Pr > -7 dB. Either one of these figures in combina- 
tion with the relative power curves will allow a determination of the relative 
frequency of fading for the rough earth case. 

3. 1.3 Surface Reflectivity Distribution 


For reflection of a plane wave by a smooth surface, most of the 
reflected signal Comes from the first few Fresnel zones. In fact, the 
total reflected signal is approximately half the signal reflected by the first 
Fresnel zone. For small transmitter and/or rooeiver heights above the 
reflecting surface, other important contributions in the reflected field come 
from the area just in front of the transmitter and receiver. For the LAS/DRS 
geometry, the major reflecting area is the first few Fresnel zones. 

As an illustration of the size of the active scattering region, 

Table 3-1 shows the size of the first Fresnel zone for various conditions. 


RELATIVE MAGNITUDES OF REFLECTIVITY DISTRIBUTION 



Figure 3-10. Reflectivity 
Distribution Near 
Transmission Plane 


Figure 3-11. Reflectivity 
Distribution Normal to 
Transmission Plane From 
Specular Point 
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TABLE 3-1. SIZE OF FIRST FRESNEL ZONE 
(LAS altitude = 300 n, mi, ) 


■ ' — ■ L 

Frequency, MHz 

U5 

400 

8000 

Satellite separation angle, 
<p in degrees 

0 

90 

0 


0 

90 

First Fresnel zone area, 
square miles 

0.4 

2. 0 




0.033 


Table 3-1 emphasizes how extremely small the active scattering 
region really is for the perfectly smooth earth. These regions correspond 
to a subtending central angle less than 0. 01 degree. 

For a rough surface, most of the reflected energy comes from near 
the specular point, but the active scattering region is much larger than the 
Fresnel zones of the smooth earth model. The surface roughness, deter- 
mines the size of these regions, whereas, in the smooth earth case, 
optical considerations determine the Fresnel zones. In addition, the 
antenna gains in the direction of a scattering point and the corresponding 
reflection coefficient also affect the power reflected from the surface near 
that point. Thus, both the gains and polarizations_for each antenna and for 
every scattering point must be known in order to calculate the relative 
power reflected from each point. 

However, to provide some insight into the reflectivity distribution, 
both the antenna characteristics and the reflection coefficient will be dis- 
regarded by defining a geo metr ical reflection factor, Fp, given by 


V 


D 


S(0) 


R 


4 ir V, 


D 


2 4 v 

r) cos Y 


exp - 


tan 2 Y ’ 


dS 


This factor is a function of the separation angle, LAS altitude, and surface 
location. If the altitude is kept constant, variation with <p and distance 
from the specular point can be shown graphically. Figure 3-10 shows F.. 
as a function of the distance from the specular point along a line parallel* 
to and very near the transmission plane, the distance from it being about 
15 miles. Figure 3-11 shows Fp as a function of the distance from the 
specular point in a direction no-mal to the transmission plane. 

Note in Figure 3-10 that the surface distribution is significantly 
different for </> - 90 degrees than for the smaller separation angles shown. 
This result is consistent with the- analytical results of Spizzichino 
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(Reference 2, pp. 257-264) showing that, for small grazing angles, a 
significant portion of the reflected energy comes from points nearer to the 
receiver and transmitter than the specular point. Small grazing angles 
resulting in this phenomenon correspond to large separation angles. 


3.2 TIME RESPONSE 

Since, as seen in Subsection 3. 1. 3, the scattering region is a signi- 
ficant geometric area on the earth's surface, energy reflected toward the 
receiver from different parts of this region arrives at different times. 
Further, from the geometric considerations of Subsection 2. 1, the length of 
the direct path and specular point reflection path are different for various 
values of separation angle. 

A simple calculation involving path length shows that,, for low-altitude 
orbits, i.e. , altitudes- between 100 and 1000 miles, the direct path time delay 
lies in the range of 115 to 155 milliseconds. Both the direct path- time delay 
and time delay for specular point reflection as a function of separation angle 
can be computed and plotted. But a more meaningful quantity for this study 
is the difference of these two time delays. This difference, defined as the 
specular differential time delay, is shown in Figure 3-12 as a function of 
separation angle, 6. The reflection path is always longer than the direct path, 
and so the specular differential time delay is the amount of time the energy 
reflected from the specular point lags the corresponding direct path energy. 

It can be shown that the specular reflection path, i. e. , VLS I v DS 
in Figure 2-4, is the shortest of all possible- earth reflection paths. Thus, 
energy reflected to the receiver from other points in. the scattering region 
will be delayed even longer than energy from the specular point. The 
numerical computation program of Appendix A calculates the power scattered 
from each surface patch and the associated time delay and then adds all the 
pbwer delayed within specified time intervals. These summations are 
normalized so that the total reflected power is unity. The resulting time 
spreading distribution represents the impulse response of the reflection 
process. 

The results are shown in Figure 3-13 for horizontal polarization, and 
the time response for vertical polarization is nearly identical. In Fig- 
ures 3-l4a and b, the curves have been normalized to the specular point 
power density, which emphasizes the relative shapes of the curves. In Fig- 
ure 3-13, the roughness factor, rj, is 0. 1, whereas in Figures 3-l4a and b, 

>1 = 0.05. Note that the roughness factor makes a considerable difference in 
the time spreading. Figure 3- 14a shows one curve corresponding to<|> = 0 
and n - 0.02, resulting in very small time spreading. For 0 - 0.05, time 
spreading is on the order of 30 microseconds, and for q = 0.1, the spread- 
ing is 80 to 100 microseconds. 
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Figure 3-14. Normalized Power Time Response 
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One of the more interesting observations concerning these figures is 
that time spreading is slightly larger for a separation angle of 30 degrees 
than for either 0 to 60 degrees. This occurs because there are two opposing 
mechanisms at work in the time spreading process: 

1) As <p increases from zero, the scattering region increases (see 
Subsection 3. 1 _3-)» tending to cause a greater time spreading. 

2) As <p increases from zero, the difference between individual 
path lengths decreases, tending to decrease the time spreading. 

Near y = 0 mechanism 1 dominates, while for large <p, mechanism 2 is 
dominant. 


3.3 FREQUENCY EFFECTS 

The frequency effects are due to the motion of the two satellites, 
which causes a doppler shift of both the direct and reflected signals. The 
frequency shifting effects may be neparated into two components as was 
done for the time response: 1) doppler shifts of the direct and specular 
point paths, and 2) frequency spreading of the reflected signal. 

The relative orientation of the two velocity vectors is the most 
significant variable- To illustrate the frequency effects and to provide 
bounds for these effects, three geometrical cases were chosen: 

Case 1: The velocity vectors of both the LAS and DRS lie in the 

transmission plane and are pointed toward each other. This 
is the most extreme case and, assuming that the DRS is in 
an equatorial orbit, would require that the LAS also have 
an equatorial orbit with opposite orbital motion. 

Case 2: The velocity vector of the LAS is in the transmission plane, 

while the DRS velocity is normal to it. This case corres- 
ponds to an LAS polar orbit and the DRS lying in the LAS 
orbit plane. 

Case 3: The velocity vectors of both satellites are normal to the 

transmission plane and lie in opposite directions. 

These three cases are illustrated in Figure 3-15, 

The change in frequency of a received single frequency signal when 
there is relative motion between the transmitter and receiver is given by 
the simple formula 


Af 


£ 

c 


dV 
o 

dt 
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Figure 3*15. Three Cases of Relative 
Velocity Orientation 



Figure 3-16. Relative Specular Doppler Shift 
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whore f is the transmitted frequency,, c is the speed of light, and V 0 is the 
distance between the two satellites. Thus, if V Q is increasing, the received 
frequency is less than the transmitted frequency. For the reflection path, 
the rate of change of the path to the earth and from the earth must be added. 
The relative doppler shift is defined as the fractional change in frequency, 
Af/f. 


Figure 3=16 shows the relative doppler shift of both the direct and 
specular reflection paths for Case 1 1 and an LAS altitude of 900 miles. Of 
particular interest in analyzing the communication problem due to multipath 
interference is the difference between these two doppler shifts,. This differ- 
ence is shown in Figure 3 = 16 and again with an expanded scale in Figure 3=17. 
Figure 3-17 also shows the relative differential doppler for Case 2 and a 
modified Case 2 corresponding to the .LAS velocity vector making a 45-degree 
angle with the transmission plane. Figure 3-18 shows the relationship 
between the relative specular differential doppler shift and altitude for 
Cast 1 2. Figure 3-19 shows the actual doppler shift in Hertz for Case 2 at 
140 MHz and also illustrates the effect of LAS altitude. 

Since the scattering region can be visualized as many small contig- 
uous scattering patches, the reflection path to a given patch will differ from 
other paths. Consequently, the doppler shift will be different for each 
patch, with a received reflected frequency spectrum being the Composite 
effect upon a single transmitted frequency.. Since for LAS/DRS geometry, 
the scattering region is nearly symmetrical about the specular point, it is 
expected that the received reflected spectrum will be nearly symmetric 
about the power density from the specular point. 

Figure 3-20 presents the reflected spectra normalized so that the 
value of the specular point is unity. This figure illustrates the relative 
doppler shift of both the direct and reflected signals, as well as the spectrum 
spreading of the reflected signal. Increasing the altitude and/or increasing 
the roughness parameter, q , would result in greater spreading of the 
spectra . 

Two points concerning Figure 3-20 are worth noting: 1) frequency 

spreading decreases as the separation angle, <p, increases, and 2) for the 
parameters indicated in the figure, the frequency of the direct signal does 
not fall within the spectrum of the reflected signal for separation angles 
greater than 30 degrees. 

Figure 3-21 illustrates the received spectra for velocity orientations 
corresponding to Case 3. For this case, the specular differential doppler 
shift is zero, i. e. , there is no doppler shift between the dirert signal and 
the signal reflected from the specular point. Thus, the only effect is the 
frequency spreading shown in Figure 3-21. So the direct signal frequency 
lies at the center of these spectra. 
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Figure 3-17. Relative Specular 
Differential Doppler Shift 
Versus Separation Angle, 
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Figure 3-22, Relative Bandwidth of 
Received Reflected Spectra 

Figure 3-22 shows the 3-dB relative bandwidth of the received 
reflected spectra as a function of the separation angle. The band' widths 
for two roughness factors are shown, and, as expected, the bandwidth from 
the smoother surface is less than that from the rougher surface. This fig- 
ure may be combined with a plot of the relative specular differential doppler 
to estimate the range of separation angle for which the relative speculat- 
differ ential doppler shift is greater than one-half the reflected spectrum 
bandwidth, i. e. , where the direct signal frequency lies outside the major 
portion of the spectrum of the reflected signal. Figure 3-23 is such a 
combination of these two curves for Case 2 and an LAS altitude of 300 n. rru. 
It can be seen that, for a roughness factor of q = 0. 05, the direct signal lies 
outside the bandwidth of the reflected spectrum for 20 < q> < 103 degrees. 

If rj =0. 02, this occurs for <p > 10 degrees. 
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Figure 3-23. Direct Signal Frequency Relationship to Reflected Spectrum 
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4. COMMUNICATION IMPROVEMENT TECHNIQUES 


In Section 3, the various aspects of the multipath phenomenon have 
been quantitatively displayed and discussed. 4 he next Logical considc l at ion 
is the manner in which communication in the multipath environment can be 
improved. The link between the LAS and DUS consists of the transmitter 
and receiver and the environment between them called the "channel. " I* or 
the sake of this discussion, the channel includes the antennas and feeds, 
while receiver and transmitter refer, respectively, to the initial and final 
electronic stages. 


From a general viewpoint, assuming that the best possible receiver is 
used, there are three basic approaches to communication improvement: 

1) Increased transmitter power 

2) Improved channel characteristics 

31 Signal processing 

The first approach has limitations which are somewhat intuitive; it 
more power is transmitted, then the inter for ring, reflected power will 
increase with the direct power. Increasing power does provide some 
improvement, but because the reflected power is also increased, this 
method is inefficient. 


The second approach refers to antenna and feed designs. The physical 
medium butvvuun two antennas is determined by the orbits of the two 
satellites, but the antennas thomsrlvrs, as part ol the channel, may he 
designed to discriminate* against the reflected signal. Both the directivity 
and polarization of the antenna may be used to improve communication. 

These techniques are discussed below, but since the emphasis in this 
study is on a low data rate omnidirectional (LAS) link, the emphasis is on 
pola rization discrimination. 


The third approach concerns the manipulation of the signal itself 
tin* dire v i signal iiiay be disv r irnmated liom the ri fleeted signal at 


so 


I t t a t I 11V VII. A V *-> *■ *- - / — 

the receiver. Almost all methods in this category result in an increased 
transmission bandwidth. Feasible techniques include pseudo-noise coding. 
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diversity techniques, frequency hopping, burst transmission, and data rate 
limiting. Plies e topics are discussed very briefly later in this section. 


4. 1 DIRECTIVE ANTENNAS 


Ihe interference due to the reflected signal can be essentially 
eliminated by using directive antennas. The purpose of this discussion is 
merely to indicate the effectiveness of this approach. Figures 4- la aftd b 
show the angular separation between the direct transmission path and the 
reflection path to the specular point as a function of the separation angle, a 
for the LAS and DRS, respectively. From Figure 4- lb, it can be seen that 
ihe DRS antenna must be highly directional in order to discriminate between 
the direct and reflection paths. Combining Figures 4- lb and 2-9, it can be 
seen that, for a 500-mile polar orbit (i = 90 degrees), the direct reflection 
path separation is greater than 1.0 degree about 70 percent of the time, but 
is less than 1. 0 degree 30 percent of the time and is always less than 
l . 8 degrees. 


A DRS antenna with a 2. 0-degree beamwidth would result in a 3-dB 
reduction of multipath interference 70 percent of the time. Although 
thi-s is helpful, it does not yield the significant advantage that a les^ 
directive antenna on the LAS will give. 

Consider a planar array with a gain of 20 dB which corresponds to 
a beamwidth of about 17 degrees. Then for a 500-mile polar orbit a 3-dB 
improvement is possible 90 percent of the time. The gain of the antenna 
drops below -15 dB for angular deviations from the boresight greater than 
1 , _ degrees. From Figures 4-la and 2-9, such an antenna will result in 
*‘r on" d ,« im P roVement 80 percent of the time. If the LAS antenna has a gain 
of 30 dB with a 5.0-degree beamwidth, a 15-dB improvement over ihe 
omnidirectional case is possible 93 percent of the time for a 500-mile 
orbit. 


The above discussion indicates the effectiveness of directive antennas 
particularly on the LAS, in reducing the interfering multipath signal. For 
any given antenna whose radiation characteristics are known, Figure 4- la, 
and 2-9 can be used to evaluate the improvement in multipath reduction. 

4. K 1 Broad Coverage Antenna 


The use of a broad coverage, but not omnidirectional, antenna may 
meet mission requirements while reducing multipath interference. An 
idealized broad coverage antenna pattern is shown in Figure 4-2 along with 
tlie resultant received relative power, R p , for horizontal polarization and 
sea ret lection at 13^ MIIz. The antenna is mounted on the LAS so that the 
maximum gain axis, X, coincides with L| in Figure 2-3 (the local vertical). 
Ihe gain pattern is a figure of revolution about the X-axis, Note that 
because of the null in the pattern in the -X - - L, direct ion, the pattern 



reduces the multipath signal the most •vhfrc 
worst — at <p - 0. 


it would otherwise be the 


4.Z POLARIZATION DISCRIMINATION 

Tn poneral the LAS spacecraft transmits an elliptically polarized 
wave which is received by the DRS with an elliptically 

Ac cases either or both spaceeralt may have lineally polarized 

antZaT. Should both spacecraft 

Reived, and, hence, 

So c P o"“atio„ would bo possible. For 

hfv?a linear antenna aid the other a circularly polarised antenna, contmua.- 
cation would be possible whatever the orientations of the antennas, and _ - 

relative power would be essentially that shown in Figure. 3-2a and b. How- 
ever, elliptically polarized antennas may be employed on bolh spacecra 
as to favor the direct signal over the reflected signal. 

Polarization discrimination is based on the fact that, during rejection, 
,he sense of rotation of a circularly polarized wave ,s reversed tf the angle 
of incidence is smaller than the pseudo-Brewster angle, Referring 
Fibres 2-16 and 2-18, it can be seen that the horizontal y polarized radia- 
tion is ref Looted with a phase shift of approximately 180 degrees with ^ 

c j. i- 0 ^he vertically polarized radiation. Thus, a circular y po arize 
^ nnal horizontal and vertical components has its sense of 

^UoT™cX,i^r^lUr M«t circularly poUrized antennas have 
at least a 20-dB rejection of circularly polarized .waves of the opposite 
se^st- therefore, L antenna adjusted to receive ihe direct signal w, 11 
reject the reflected signal. 

To more clearly illustrate the effect of antenna polarization, con- 
sider the polarization coefficient, PC, given by Equation 58 in 
Subsection 2. 6. 
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where the polarization factors are defined at the beginning of Subsection 2.6. 
U the 1JUS antenna is circularly polarized, then from liquation 41 in 
Subsection 2* S 
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Similarly, il (In* LAS transmits circular polarization in both tin* direct 
path and scattering region directions, then 


For angles less than the pseudo-Brewster angle, the phase of the 
polarization is reversed, resulting in P,. j. Then, clearly, 


vertical 


t l P P 

r rr 


1 l j 2 - 0 


and so PC 


0 . 


Of course, a real LAS antenna will not transmit exact circular 
polarization in two different directions, and the polarization is not completely 
reversed during reflection. This can be accounted for by antenna ellipticity. 

The relative power expressions given by Equations 63 and 64 in 
Subsection 2.6 contain the polarization coefficient as a factor and may be 
used to evaluate the effects of elliptical polarization. It was shown in 
Subsection 3. 1 that the smooth earth model gives results very close to 
those of the rough earth model. Thus, most of the numerical computation 
associated with this section used this model because of the reduced machine 
compulation expense. Several rough earth computations were made to 

result aPpr ° aCh * FiBuri ‘ B 4 ' 3 lhrou « h Present the computation 

I* igure 4-3 shows the effect of polarization discrimination on the 
relative power as a function of separation angle for a 140-MIIz transmission 
frequency and sea water reflection. Since the direction of the scattering 
legion and the direct path is nearly the same for the DRS because of its 
relatively large dis'unce from the earth, the polarization can be eonsidered 
constant for both the direct and reflection paths. Thus, Figure 4-3a 
corresponds to a circularly polarized DRS antenna with respect to both the 
et lection and direct paths (polarization ellipticity = 0). Figure 4-3b corres- 
ponds to a DRS polarization ellipticity of -4 dB, which represents a greater 
c mount of veitieal polarization. The nine curves in each of these figures 
correspond to polarization ellipticities of the LAS antenna. Positive ellio- 
ticities repres- nt more horizontal polarization than vertical and negative 
t* 111 pfcic i tit* 8 imply the converse . 

i\o(e th.nl in 1 igure 4-3a the rvo corresponding to a circularly 
polarized LAS antenna results in excellent reduction of R„. m Figure 4-31, 
tile most Sigmlieanl reduction occurs when the FAS has the opposite type of’ 
polarization. That is, in Figure 4- 3b, the DRS is mostly vertically polarize.!; 
Huis, R the DAS antenna is mostly horizontally polarized, the relative power 
will he greatly reduced. 
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npc: n„i. Fl ^ UrC 4 -~^ V s t ^ l ° l ' osult of rou Sh earth computations for both IAS and 

For r r h ^T^;^r io; paths * 

mcnt between these results and the circular pi^^Zouot a « r ° e ’ 

rouehMHh 8 ' 4 -4a, and 4-5a. For smaller separation angles, the 
lough eaith calculations are deemed to be the more correct. 

4. 2. 1 Conclusion 


r rom examination of Fiouro 4-S smrl a r ^ 4l 

mission frequencv of 140 1V114V ’ j i i ..." *, a PP ea rs that, at a trans- 

.he LAS It^theT,^ “ 
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SIGNAL PROCESSING 


E»F#i'=5S“^S£ S&r 

4 « 3 * 1 Pseudo-Noise C oding 

called maxVnanO^gVh^ ° f PN se< <" e,u ' , ‘ 3 

1) In each period of the sequence, the number of ONE's differs 

from the number of ZERO'S by, at most, one. ‘ 

2) Among the runs of ONE's ana of ZERO's in each period, one- 
huff ol each kind are oi length one, one-fourth of each 
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Figure 4-5. Polarization Discrimination for Smooth Earth Model 
Land reflection for 140 to 2 GHz 
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Figure 4-6. Polarization Discrimination for Rough Earth Model 
LAS and DRS antennas circularly polarized 
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kind are of length two, one-eighth are of length three, and ko 
on as long as these fractions give meaningful numbers of runs, 


3) Ii a period of the sequence is added modulo 2 with any cyclic 
shift itself, the number of agreements differs from the number 
of disagreements by one at most. 


This last propet ty can be restated as follows: For integral phase 

displacements, the autocorrelation furction of a PN sequence is two valued, 
Specifically, it looks as is shown in Figure 4-7, where, if s(t) represents 
a normalized PN code signal 


.T 


r 6 m - T 


S*(l) s(t I T) dt 


It is the nature of the autocorrelation function which makes a PN sequence 
useful for signal detection. 

13y performing the correlation process on the received signal and 
detecting the peak of the autocorrelation function, interfering noise and 
signals such as the multipath signal can be rejected. One requirement for 
a PN system to operate properly in the LAS/DUS multipath environment is 
that a code pulse length be shorter than the delay between the direct and 
reflected signals. Referring to Figure 3-13, if the code pulse width is 
0. 1 millisecond, then this system will operate correctly in a 300-mile 
altitude satellite for separation angles less than 93 degrees. When the sepa 
ration angle is larger, the time delay between the direct and reflected 
signals is Less than the pulse width which corresponds to an RF bandwidth 
ol approximately 30 kHz. Decreasing the pulse width, which requires an 
increase in bandwidth, will allow preater communication time. Thus, 
mission requirements and frequency allocations vv : ll have a major influence 
in determining the code pulse width. 

4. 3. 2 Diversity Techniq ues 

Diversity techniques involve the establishment of N distinguishable, 
dissimilarly fading signal transmission channels. The diversity receiving 
system chooses at each instant a desirable combination of these signals. 

The techniques include space diversity, angle of arrival diversity, polar- 
ization diversity, frequency diversity, time diversity, and multipath 
diversity, lie* >irst three are not purely signal processing techniques as 
defined above. Space diversity and angle of arrival diversity require two ot 
more antennas and receivers, which is somewhat impractical on a spacer ra 
Polarization diversity requires crossed linearly polarized feeds wifh 
separate receivers. Frequency diversity requires transmission of the sign; 















on two or more different carrier or subcarricr frequencies. Time diversity 
and multipath diversity are closely related and have been applied only to 
digital transmission. In time diversity, the same information bit is 
repeated at time intervals greater than the reciprocal of the fading rate. 

This technique requires storage at both the transmitter and receiver. If 
the multipath consists of many channels with various delay times, with a 
multipath diversity system the received signal can be processed to separate 
the signals corresponding to different paths and then be recombined. These 
techniques, along with combination analysis are treated at length in 
Reference 24, 

Diversity techniques are usually employed in a multipath environ- 
ment similar to the reflection path alone, and they are used to combat the 
Rayleigh fading that occurs in this phenomenon. But, for the geometry and 
physics of this problem, there are only two basic paths rather than many, 
even though the reflection path consists of a collection of many subpaths. 

Thus, a diversity technique would establish a number of channels to provide 
a signal combination which would ameliorate the fading caused by interfer- 
ence of the two basic transmission paths . 

4. 3. 3 Frequency Hopping 

In the frequency hopping scheme, the carrier frequency of the trans- 
mitter is switched in a cyclic progression through several values. The 
receiver must be synchronized to the transmitter and thus rejects the delayed 
multipath signal whose frequency differs from that being instantaneously 
accepted by the receiver. To clarify this technique, consider an example 
where each frequency is to be transmitted for a period equal to one-fifth 
the maximum specular differential time delay, T c { m (see Subsection 3.2). 
After completion of transmission at one frequency, the corresponding multi- 
path signal may arrive at the receiver until T clm seconds after completion 
of the transmission. Thus, 1. 2 Tq m seconds are required between the 
beginnings of transmissions at each frequency, and, since each transmission 
lasts 0.2 T r j m seconds, six frequencies will be required. Figure 4-8 
illustrates these concepts. 

Power requirements are the same as for continuous transmission at 
a single frequency, but the bandwidth must be increased to something larger 
than six times the original bandwidth. The relationship between the original 
bandwidth and the frequency hopping bandwidth depends on the size of the 
original bandwidth and the doppler shift effects (see Subsection 3.3). 

It should be noted that, in the above example, multipath inter- 
ference can be eliminated when the multipath delay is greater than 0. 2 T . , 

but, when it is less, interference will occur. If the LAS altitude is 500 * m 
miles, then according to Figure 3-12, the system in the above example 
would not eliminate interference for separation angles larger than about 
80 degrees. Thus, from Figure 2-9, this system would eliminate inter- 
ference about 54 percent of the mission lifetime for an LAS polar orbit, and 
about 74 percent of the time for an equatorial orbit. 
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4. 3. 4 Buret Transmission 


Rotated to the frequency hopping method is the technique of burst 
transmission. Only one transmission frequency is used, but data is trans" 
mitted in periodic bursts where a burst transmission period is a fraction a 
of the maximum specular differential time delay, T c j m . The time between 
bursts is (1 f «)T f i m , and so the data compression factor is given by 
G t i ~ (1 +a)/ft. If a is again chosen to boQ.2, then C ( j s 6, and the RF 
bandwidth must be increased by this factor. For this system, the peak 
power increases also by the factor C f j, but the average power is the same 
as for continuous transmission. 

4. 3. 5 Data Rate Limiting 


With data rate limiting, a simple technique, data bits are transmitted 
with pulse lengths which are large compared to the maximum specular dif- 
ferential time delay, T c i m . Thus, if the pulses have lengths of PT^ni* where 
(3 is a factor based on system considerations ((3 > 2), then the bit rate, R^ is 
given by 


R 


b 


1 

P ^dm 


If p s 4 and the LAS altitude is 500 miles, then from Figure 3-12, 

= 46. 5 bits/ sec. This rate may be adequate for housekeeping telemetry 
and command data. For a given value of /3, the data rate will decrease 
with increasing LAS altitude. 


5. EXPERIMENTS 


Experimental data is needed to validate the statistical and 
mathematical assumptions employed in developing the mathematical model 
of Subsection 2. 4. Sections 3 and 4 present the results of computation 
based on this model, but this quantitative data is no more valid than the 
assumptions used to produce it. If experimental data differs significantly 
from the results based on the model, then this data may be used to modify 
the model, i. e. , to correct and calibrate it. Whether or not the experi- 
mental data agrees with mathematical results, it is required for.oon£idence 
in the characteristics of the multipath phenomenon. 


5. 1 EXPERIMENT-DEFINITION 

The characteristics to be measured experimentally are those dis- 
cussed in Section 3, namely the relative power, time response, and fre- 
quency effects. Of these three, the relative power, i.e., the ratio of 
average reflected power to direct path power, is most important. The time 
spreading of the reflected signal is much smaller than the differential delay, 
which can be computed from purely geometric considerations, and,, hence, 
is of least importance. Measurement of frequency spreading will allow 
evaluatio i of the statistical assumptions used in developing the rough earth 
mathematical model (see Subsection 2.4). 


A number of propagation aspects influence the multipath phenomenon; 
they are quantitatively included in the mathematical model- in the parameters 
of the basic formulas (see Section 2). These aspects include: 


Surface state: Determines the roughness factor, n, of Equations 63, 

77, and 78 


Surface type : The electrical properties of the surface influence the 

reflection coefficient, R c , in Equations 20 and 37 

Pola rization : The polarization of the electromagnetic wave with 

respect to the surface and the receiving antenna 
influences the reflection coefficient and the polariza- 
tion efficiency, respectively (see Subsections 2.4 
and 2. 5) 
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I ransmission Intiuencps both t lie reflection coefficient and the 
frequency ; smoothness factor 

Geometry ; This is probably the most important aspect since 

the separation angle i,s the single most important 
parameter. The altitudes of the satellites also 
influence the reflection and direct path character- 
istics. In particular, the distances V u , \Cj v and 
Vf)! the reflection coefficient; and the incidence 
angle, 0, are a function of the geometry. For known 
antenna characteristics, l'*Uh the gain and polariza- 
tion for the two paths will be determined by the 
geometry. 

A complete experimental program will provide a sufficient variation 
of each parameter associated with these aspects to evaluate the validity of 
the mathematical model over the range ol interest. For the roughness 
factor, r|, a range of variation should be established. 

5.1.1 Signal Reception 


The evaluation of data is facilitated if signals from tile reflection 
path and direct path can be received separately. This can be accomplished 
by either two antenna and receiving systems or a single receiving system 
which separates reception from the two paths by time. That is, first the 
reflected signal is passed through the receiver and then processed, and then 
the direct signal is passed through the same receiver and processor. The 
scheme employing two receivers eliminates any possibility of error due to 
the time separation of the two measurements, while the latter method, 
employing only one receiver, eliminates the need for calibration. 

This separation of signals requires directive antennas. From Fig- 
ure 4-1, itcanbeseen that,, for a satellite -to -satellite link, beamwidths as 
large as 20 degrees will allow signal separation over a wide range of separ- 
ation angles. However, at aircraft altitudes, the angular separation of the 
two paths is much smaller, and, consequently, a highly directive antenna 
will be required to separate the two signals. Directive antennas imply the 
use of higher frequencies, yet some of the multipath characteristics are 
particularly of interest at VHF where directive antennas arc not feasible. 
Thus, while many aspects of the mathematical model may be checked with 
the use of high-frequency transmission (I GHz and above), a complete evalu- 
ation ol tilt' model will require measurements at VHF with broad coverage 
antennas. 

By proper shielding, VHF antennas can be made somewhat directive, 
thus effecting a separation of the reflected and direct signals over a range of 
geometric variation. But when separation of the two signals with directive 
antennas is not leasiblc or practical, modulation and signal processing may 
be employed to aid in evaluating the multipath characteristics. Further, even 
with separated signals, the modulation technique affects the ease ol' process- 
ing and evaluating the data. 
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5. 1 . Modulation 


An attempt wan made to evaluate the- efficiency ot various types of 
modulation in determining various cha raete ri i sties of the reflected signal 
when the signals ran ho separated* These qualitative results are shown in 
Table 5-1 , Thun, when rotTeeted and direct signals ran be separated, trans- 
mlHHion of a single frequeuey rarrier will allow determination oJ‘ the major 
eha raete ri sties, i. e. , the frequency spreading, polarization elTe ets, reflec- 
tion coefficients, and hence, the relative power. The time delay and time 
sp reading may be measured by short pulses* 


When the two signals cannot be* separated, short pulses again will 
allow evaluation of all characteristics., but the processing complexity is 
much greater than for the separate signal case. 


5. I LABORATORY MEASUREMENTS 

A relatively simple and inexpensive series of experiments could be 
performed using simulated sea surfaces in a laboratory environment. Based 
upon existing knowledge of sea sta.o and si' a statistics, a computer -gene rated 
sea surface may be constructed of clay or plaster of paris and then given the 
appropriate electrical properties by spray-coating this surface. This type of 
surface model would allow variation of the geometry, polarization, and trans- 
mission frequency, and frequency spreading and relative power could be 
evaluated. These two reflection characteristics may be used to evaluate the 
validity of many of the assumptions and approximations used in the develop- 
ment of the mathematical model. Further, several surfaces could be simu- 
lated in order to allow variation of the roughness and electrical properties. 


TABLE 5-1. TYPE OF MODULATION* 



Data That Can Ik- Extracted From Multipath Signal 

Type of 
Modulation 

Differential 
Time Delay 

Frequency 

Spread 

Polarization 

Reflection 
Cha raete risties 

cw 

No 

Yes 1 

Yes 1 

Yes 1 . 

Narrow pulses 

Yes i. 

Y e s 1 

Yen L 

Yes t! 

FM 

Yes 4 

Y e s 3 

Yen 3 

Yes 3 

Frequency 

sweep 

Yes 1 

Yes 3 

Yes 1 

Yu s 3 

PCM - FM 

Yes 1 

i 

Yes 3 

Yes 3 

Yes 3 


•"The numbers 1 through 4 indicate the range of implementation difficulty, 
i. r. , 1 easy, 4 difficult. 
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REAL ENVIRONMENT EXPERIMENTS 


5. 3 


Although laboratory measurements can provide an intermediate and 
useful step in evaluating the charade risties of the multipath phenomenon and 
in evaluating the mathematical model developed in this study, confidence 
will he lacking in the conclusions drawn from these results because the mea- 
surements are based upon a simulation of the natural phenomenon rather 
than the actual phenomenon itself. It is, therefore, most natural to turn to 
the real earth reflection process in order to gain confidence in the data and 
the mathematical model it may validate. Feasible experiments, in order of 
sophistication and usefulness, are described below. 

5.3.1 Ai rcraft/Aircraft Propagation 

Transmission between two aircraft with highly directive antennas 
will allow reflection coefficients and polarization effects to be evaluated. 
Relative power and frequency spreading will provide data points for verify- 
ing the mathematical model. The geometry, frequency, and polarization may 
be varied, and it is relatively easy to make measurements of wive height in 
the scattering region. 

5.3.3 ATS- V/Aire raft 


A directional antenna on an aircraft, which may bo pointed upward to 
receive the pulse -like transmission direct from the ATS-V spinning space- 
craft and then toward the earth for receiving the reflected signal, Would 
allow determination of relative power for verifying the mathematical model. 
Because of the short pulse nature of the received transmission, frequency 
spreading data will require more processing than the aircraft/aircraft 
experiment. The advantages of this experiment over that described in Sub- 
section 5.3. ] are: 1 ) only one aircraft need be equipped and flown, and 

2) the geometry is closer to that of the LAS/DRSS geometry, resulting in 
scattering from a 1 >er region than the aircraft/aireraft experiment. The 
disadvantages arc: i, .he frequency is fixed in the L-band, and 2) polariza- 
tion measurements are limited. This experiment is particularly attractive 
because of the relative ease of implementation and the near LAS/DRS geometry. 

5.3.3 ATS-F/Nimbus E 

At present, an experiment is planned for data transmission between 
the Nimbus E, low-altitude spacecraft, and the ATS-F synchronous satellite. 
Thus, the geometry is exact 1 y that of interest in this study. An S-band direc- 
tive antenna is used for this link. The feasibility of using these two satellites 
for multipath measurements depends on the ability of the Nimbus E to point 
the directive antenna toward the earth. At the time of this writing, the 
antenna's location and pointing capabilities allow the antenna to point at the 
earth only for large separation angles when Nimbus E is operating in its 
normal earth pointing mode. If the attitude could be changed and monitored, 
the antenna could point toward the specular point for earth reflection to 
ATS-F. This experiment would provide excellent data and allow an evalua- 
tion of the actual effects experienced in the LAS/DRSS environment. 
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5.3.4 ATS-G/LAS 


An experiment specifically for measuring multipath characteristics 
could bo designed to-employ a low -altitude satellite transmitting to the 
ATS-G spacecraft. This experiment should be designed to supplement 
data acquit ed by any previous experiments and, further, should 
allow variation of as many parameters as possible while providing data for 
evaluating the three effects mentioned in Subsection 5, 1, 
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APPENDIX A. COMPUTER PROGRAM FOR 
MULT IPA Til CO MP UT AT 10 NS 


A computer program was written in the FORTRAN IV language for 
computing the quantities of interest in Section 3, i.e,., total relative power, 
time spreading, and frequency spreading. The basis for computations is 
given in Section 2, and the final equations which are employed in the pro- 
gram are given by Equations 63, 64, 77, 78, 80, and 81 of Subsection 2.6. 
There are many quantities in the expressions of these equations which must 
be determined in intermediate steps. 

The only geometric input quantity is the LAS altitude. From this, 
the maximum value of the separation angle, <)>, is determined, and then the 
program progresses in increments of 4> from an initial value, both of which 
arc input quantities, to a final value of <t>, which differs from<j) max by less 
than the increment. For each value of <}>, a complete computation of all the 
characteristics is performed. 

The specular point is determined first by an iterative procedure, and 
then many quantities related to the specular point are calculated. These 
include : 

1) Reflection coefficients — horizontal-and vert 4 al 

2) Geometric angles shown in Figure A-l 

3) The magnitudes of Vj^, and Vpj shown in Figure A-l 

4) Specular time delay 

5) Differential time delay 

6) Specular relative dapple r 

7) Differential relative doppler 

8) The differential doppler for the input frequency of interest 

V) Coherency factor 
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IQ) Divergence factor 
J J ) Shadowing coefficient 

The input quantities whieh are used in the above computations are: 

1) Surface relative dielectric constant 

2) Surface conductivity 

3) Transmission frequency 

4) Roughness factor, q 

5) RMS surface variation, tr 

Characteristics for antennas different from the two linear orientations 
defined in Subsection 2,6 can be specified via subroutines. The output 
of the subroutines are gains and polarization factors in the direct and scat- 
tering point directions for both the DRS and LAS antennas. -Using these, the 
program also computes a value for the polarization coefficient at the- specular 
point. 


With the above computed quantities, the relative power is computed 
for a perfectly smooth earth for the two linear polarizations and for the sub- 
routine antenna characteristics. Using the coherency factor corresponding 
to the input frequency and rms surface variation, the smooth earth relative 
power multiplied by the coherency factor yields an estimate of the coherent 
relative power. The program then begins its rough earth calculations. 

The rough earth relative power estimates are a result of numerical 
integration over the earth' s surface. However,, since the significant reflected 
power is scattered from a limited region surrounding the specular point, -it 
is ncCessai. y to determine some bounds on. this region. Referring to Fig- 
ure 2-7, the integration is performed for surface increments dS by varying 
the two geocentric angles a and p in small increments. The size of a patch, 
is RE^ cos p do dp. The size of the scattering region varies with the separa- 
tion angle as shoWn in Subsection 3. 1. 3. But, limits on the angles a and p 
can be established for a given value of <p. It can be shown that for <p - 0 


V 


ts 


R„ d a 
2h 


where h represents the LAS attitude. 

Using the above equation, define a a by 

2 h n 

Q —s — L 

° r e 
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which in the value of p where v «'|* Since the scattering crons section is 
a function of f'stpt-.tan^Y/'l^) e espi "Y'Vq*)* the Umit on both p and p for 
^ ' 0 was chosen to bo 4p (J . Hut, for 40 s <p < SO decrees, the limit on p 
was increased to 0 o u , and for HO *<p < 120 degrees, PHn>. 12 p u . The 
limit on p remained the same for all values of (p. 

To determine the integration increments dp and dp, several values 
were tried, and reliable results with reasonable computation time were 
achieved with 


d« dp 0.4 a 

o 

1 he program then pertormed the summation u£ incremental relative power 
for 


"f’lini £ P £ Plim H 4 % 

with o- and P varied in increments of 0.4 e u . For <p 0, this amounts to 
dividing the scattering region into 400 scattering patches, and for 
2 80 degrees, 1200 patches comprise the scattering region. 

The formulas used are given by Equations 63, 77, and 78 of Subsec- 
tion 2.6 where the following common factor is first computed. 


cy®, p) 



S(0) 


2 4“ e *P 

tit q cos q 



Then the factors accounting for the antenna characteristics are employed, 
tor instance, from Equation 77 of Subsection 2,6 for horizontal polarization 

(UR p’[!ll | e H1U1 (o ’ P'C P» 

anti hu 


(R p , im EKnui-P)|^s- P ) 
s 


whoro Ih glvtni by Equation 73 of that name* eubunutum. Similarly, for 

vortical polarisation, a a clofimul In Subamtlon 2.6 


(<u Vvv 


’VliVl t . 

|!»H b V 


For a general act of antenna eharaeteriatieH, 
tion 2. 6 


from Equation 63 of Subset; = 


dll 


Cl 3 G 4 


P G 1 C, > 


K, . ( P. K , 

hh i vh 


(PC) C 


where the quantities in the above' expression are defined in Section 2. 

For each value of & and p, the differential time del: V tween tie- 

reflection paths corresponding to specular point an 1 tP Scattering paten 
center is calculated. The incremental relative puv. o r is then accumulated in 
time slots, resulting in an impulse response for the reflection process. 

Similarly, for relative velocities of the two satellites, corresponding 
to Case 2 of Subsection 3,3, the relative differential doppler shift between 
the specular point uid scattering patch center is calculated. Then the incre- 
mental relative power, is accumulated in relative frequency slots, resulting 
in a relative power spect rum corresponding to a single, constant, transmitted 
frequency. However, since relative frequency is merely the fractional cliang 
in frequency, the results can be used tu estimate the effect of reflection on a 
transmitted frequency spectrum. 

The results of total relative power, time response, and frequency 
response, as calculated by thjjs^program, are- discussed in Section 3. 
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